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ABSTRACT 

We present the results of a new spectroscopic study of Fe K-band absorption in Active 
Galactic Nuclei (AGN). Using data obtained from the Suzaku public archive we have 
performed a statistically driven blind search for Fexxv Hea and/or Fexxvi Lya 
absorption lines in a large sample of 51 type 1.0 — 1.9 AGN. Through extensive 
Monte Carlo simulations we find that statistically significant absorption is detected 
at E > 6.7 keV in 20/51 sources at the Pmc ^ 95% level, which corresponds to 
<~ 40% of the total sample. In all cases, individual absorption lines are detected in- 
dependently and simultaneously amongst the two (or three) available XIS detectors 
which confirms the robustness of the line detections. The most frequently observed 
outflow phenomenology consists of two discrete absorption troughs corresponding to 
Fexxv Hea and Fexxvi Lya at a common velocity shift. From xstar fitting the 
mean column density and ionisation parameter for the FeK absorption components 
are log(iVH/cm~ 2 ) w 23 and log^/ergcms^ 1 ) w 4.5, respectively. Measured outflow 
velocities span a continuous range from < 1, 500 kms -1 up to ~ 100, 000 km s -1 , with 
mean and median values of <~ 0.1c and ~ 0.056 c, respectively. The results of this 
work are consistent with those recently obtained using XMM-Newton and indepen- 
dently provides strong evidence for the existence of very highly-ionised circumnuclear 
material in a significant fraction of both radio-quiet and radio-loud AGN in the local 
universe. 
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1 INTRODUCTION 

Observational evidence for outflows and winds in Active 
Galactic Nuclei (AGN) is seen in multiple energy regimes; 
ranging from the prominent radio-jets seen in radio-loud 
sources, to the broad absorption lines observed in Broad Ab- 
sorption Line (BAL) quasars, through to the photoionised 
'warm absorber' which is frequently observed in the soft X- 
rays (e.g., Blustin et al. 2005; McKernan et al. 2007; Cren- 
shaw et al. 2003). Gravitational micro-lensing studies have 
shown that the primary X-ray emission region in AGN is on 
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the order of a few tens of gravitational radii (R s = GM/c 2 ) 
in size (Morgan et al. 2008; Chartas et al. 2009; Dai et al. 
2010) and so the spectral analysis of absorption features im- 
printed on the X-ray continuum by circumnuclear material 
in an AGN is a powerful diagnostic of the physical conditions 
of the environment in the vicinity of the central Super Mas- 
sive Black Hole (SMBH), of the dynamics and kinematics 
of the outflowing material, its chemical composition and its 
ionisation state. Understanding how such winds are formed, 
their physical characteristics and their energetic output is of 
vital importance when it comes to studying how the inter- 
play between the accretion- and ejection-flows at small radii 
can affect the host galaxy on larger scales, (e.g., Ferrarese 
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& Merritt 2000; King & Pounds 2003; Sazonov et al. 2005; 
Ciotti et al. 2009) 

In the X-ray regime the presence of photoionised mate- 
rial in AGN is well established, with at least 50% of objects 
showing direct spectroscopic evidence for discrete absorption 
lines and photoelectric edges in their soft-band (E < 3 keV) 
spectra (e.g., Reynolds 1997; Crenshaw et al. 2003; Blustin 
et al. 2005), with typical line of sight velocities ranging 
from a few hundred to around a thousand kms -1 (Blustin 
et al. 2005; McKernan et al. 2007). Detailed modelling at 
high spectral resolution with photoionisation codes such as 
xstar (Kallman et al. 2004) often finds the soft-band ab- 
sorber to be described by column densities and ionisation 
parameters 1 in the range of log(iVH/cm~ 2 ) ~ 20 — 23 and 
log^/crgcms -1 ) ~ 0—3, respectively. These parameters im- 
ply that the warm absorbers are typically located on fairly 
large distances from the central black hole, and perhaps as- 
sociated with a wind originating from the putative parsec 
scale torus (Blustin et al. 2005) or with the latter stages of 
an accretion disc wind which has propagated out to larger 
radii (Proga & Kallman 2004) . By virtue of their low outflow 
velocities the soft X-ray warm absorber is thought to only 
have a weak feedback effect in their host galaxy. Indeed, the 
mechanical power imparted by individual warm absorption 
components is very low, typically < 0.01% of an AGN's bolo- 
metric luminosity (Lbol) (e.g., Blustin et al. 2005), which is 
significantly lower than the ~ 0.5% of L^oi thought nec- 
essary for feedback to affect the host galaxy (Hopkins & 
Elvis 2010). However, Crenshaw & Kraemer (2012) have re- 
cently shown that this ~ 0.5% threshold can be exceeded 
provided the mechanical power is integrated over all UV 
and X-ray absorption components, at least in the case of a 
few moderate-luminosity local AGN. 

More recently the higher throughput and larger effective 
area offered by XMM-Newton and Suzaku at higher X-ray 
energies (E ~ 5 — 10 keV) has shown that absorption, specif- 
ically in the form of very highly-ionised resonant absorp- 
tion lines associated with the K-shell (Is— >2p) transitions 
of Fexxv and Fexxvi, is also manifested in the hard X-ray 
spectrum of a significant fraction of local AGN. While evi- 
dence for such absorption lines was initially confined to de- 
tailed studies of individual sources (e.g., Reeves et al. 2003; 
Pounds et al. 2003; Risaliti et al. 2005; Turner et al. 2008; 
Cappi et al. 2009; Giustini et al. 2011) the recent systematic 
archival XMM-Newton study conducted by Tombesi et al. 
(2010, 2011, 2012) has shown that Fexxv- xxvi absorption 
lines are present in the X-ray spectra of > 40% of radio- 
quiet AGN in the local universe (z < 0.1). Moreover, such 
outflows have also been detected in a small sample of local 
Broad Line Radio Galaxies (BLRGs; Tombesi et al. 2010, 
hereafter 'T10B') which thus suggests that they may rep- 
resent an important addition to the commonly held AGN 
unification model (e.g., Antonucci 1993; Urry & Padovani 
1995. 

In comparison to the soft-band absorbers these hard 
X-ray absorbers generally have much more extreme param- 

The ionisation parameter is defined as £ = ^ion /nR 2 (Tarter 
et al. 1969), where L; on is the 1 — 1000 Rydberg ionising luminos- 
ity, n is the electron density and R is the distance of the ionising 
source from the absorbing clouds. 



eters, with log(./VH/cm~ 2 ) ~ 23 — 24 and log(^/ergcms _1 ) ~ 
3 — 6, and their outflow velocities relative to the host galaxy 
can reach mildly relativistic values. While alternative ex- 
planations have been posited in the literature, e.g., reso- 
nant absorption by highly-ionised material in an co-rotating 
and optically-thin plasma above the accretion disc (Gallo & 
Fabian 2011), the large inferred velocities - combined with 
the short time-scale variability sometimes exhibited by the 
absorption features - point to an origin more likely associ- 
ated with a wind which is launched from the surface of the 
accretion disc itself (e.g., Pounds et al. 2003; Reeves et al. 
2009; Gofford et al. 2011; Tombesi et al. 2012). In this sce- 
nario the inferred mass outflow rates for disc- winds are often 
comparable to those of the matter which accretes onto the 
central black hole and the consequent mechanical power can 
also be a sizeable fraction (i.e., ^ few%) of an AGN's bolo- 
metric luminosity (e.g., Chartas et al. 2002; Pounds et al. 
2003; Gibson et al. 2005; Reeves et al. 2009; Gofford et al. 
2011; Tombesi et al. 2012) making them a promising means 
of linking the small- and large-scale processes at play in 
galactic evolution over cosmic time. 

1.1 Why Suzaku? 

In this work we use archival Suzaku observations of a large 
sample of AGN to characterise the properties of highly- 
ionised outflows in the FeK band. To date, the most com- 
prehensive study of these outflows has been conducted by 
Tombesi et al. (2010, hereafter T10A) who performed a sys- 
tematic narrow-band (i.e., 3.5 — 10.5 keV) analysis of 42 
sources (with 101 observations) obtained from the XMM- 
Newton archive; using a simple baseline model to describe 
the AGN continuum in the FeK band. The baseline model 
consisted of a power-law, narrow Gaussians and, where re- 
quired, neutral absorption to approximate for any spectral 
curvature. This phenomenological approach yielded a fit to 
the 4—10 keV energy band of most sources which was suf- 
ficient to enable the systematic search of Fexxv Hea and 
Fexxvi Lya absorption lines without needing to take ad- 
ditional spectral complexity into account. Moreover, T10A 
found that this approach resulted in continuum parameters 
which were largely consistent with those found by authors 
who conducted a more thorough fit using the entire XMM- 
Newton bandpass. However, while this approach is suitable 
for most sources it is important to note that in those which 
have more complex X-ray spectra, e.g., those which are very 
heavily absorbed or those with strong hard excesses, using 
a narrow-band fit can yield a model which is a poor rep- 
resentation of the data when extrapolated to consider the 
higher and lower energies. The only way to overcome this 
limitation is by performing a detailed broadband analysis. 
Suzaku is currently the only X-ray observatory which offers 
a sufficiently broad bandpass (i.e., 0.6 — 50 keV) to allow for 
the effects of soft band absorption, the soft-excess and the 
Compton reflection component to be constrained simulta- 
neously. This makes it the ideal instrument to confidently 
constrain the underlying continuum and to robustly assess 
for the presence of highly-ionised outflows in AGN. 
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2 SAMPLE SELECTION 

The sample was selected from the Data Archive and Trans- 
mission System 2 (DARTS; Tamura et al. 2004) which con- 
tains all publicly available Suzaku observations categorised 
by the classification of the target source. An initial sample 
of AGN was drawn by positionally cross-correlating the tar- 
geted pointing coordinates of all publicly available (as of the 
end of December 2011) Suzaku observations of extragalac- 
tic compact sources against the known positions of sources 
contained in the VERONCAT catalogue of Quasars & AGN 
(Veron-Cetty & Veron 2010). The VERONCAT has an ex- 
tensive list of local AGN, all of which have been quantita- 
tively classified based upon their optical properties using the 
criteria introduced by Winkler (1992). Observations of non- 
AGN which were matched by virtue of their having a similar 
position to a known AGN on the sky, such as those of extra- 
galactic Ultra-luminous X-ray Sources (ULXs; Feng & Soria 
2011) or X-ray bright supernovae, were systematically ex- 
cluded. Only those observations with exposures long enough 
to ensure that the net source counts between 2 — lOkeV 
in the source rest-frame exceeded ~15,000 were retained 
(typical exposures were > 50 ks) such that a narrow (i.e., 
EW = 30 eV) absorption features was detectable at 95% 
from Monte Carlo simulations at source rest-frame energies 
of up to 8 — 9 keV (see Section 4.4). In the case of the high 
red-shift quasar APM 08279+5255, which shows evidence for 
absorption lines at rest-frame energies of E > lOkeV (e.g., 
Saez et al. 2009; Chartas et al. 2009; Saez & Chartas 2011) 
the measure of net source counts was instead taken for the 
entire XIS bandpass (i.e., 0.6—10 keV in the observer frame) 
because the FeK features are shifted to lower energies due 
to the high red-shift of the AGN. 

As highly-ionised outflows are thought to originate at 
relatively small distances from the central nucleus (e.g., Gof- 
ford et al. 2011; Tombesi et al. 2011, 2012, hereafter 'Til' 
and 'T12', respectively) it is important that the primary 
continuum emission from the central nucleus, rather than 
that which is reprocessed/scattered by circumnuclear ma- 
terial out of line of sight, is directly observed so that such 
outflows can be detected. To this end we exclude all Type- 
2 sources to make sure that all sources were optically-thin 
to X-rays below lOkeV. Therefore only those sources with 
a Type 1.0 — 1.9 orientation to the line of sight, as per 
the classifications listed on the NASA/IPAC Extragalactic 
Database (hereafter NED), in the VERONCAT catalogue it- 
self and through literary sources were included. Note that we 
conservatively include the radio-quiet AGN ESO 103-G035 
(z — 0.01329) which, despite being classified as Type-2.0 in 
the VERONCAT and on NED, is often regarded as a Type- 
1.9 Seyfert in the literature (e.g., Warwick et al. 1988) by 
virtue of the presence of a moderately broad Ha line in its 
optical spectrum (Phillips et al. 1979). Observation details 
for all sources included in the heterogeneously selected sam- 
ple are listed in Appendix A. There are 99 observations of 
51 AGN spanning a wide range of spectral types and radio- 
properties. As shown in Table 1 the sample is dominated by 
low-moderate inclination Seyfert 1.0 — 1.5 galaxies (28/51; 
34/51 if Narrow Line Seyfert Is are included), and contains 
comparatively few high-inclination (Type 1.8— 1.9) Seyferts 

2 http:/ /darts.jaxa.jp/astro/suzaku/ 
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Figure 1. Histogram showing the distribution of sources included 
in this work. The high red-shift QSO APM 08279+5255 (z = 3.91) 
has been omitted for scaling purposes 
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Figure 2. Histogram showing the logarithm of the total front- 
illuminated XIS 2—10 keV counts in all fitted spectra. In stacked 
spectra (sec table Al) the total co-added counts are considered 
here, rather than the counts in each individual sequence. 



(6/51). There are 6 radio-loud sources in the sample, includ- 
ing all 5 of the Broad Line Radio Galaxies (BLRGs) included 
in the T10B outflow case study, and 5 QSOs. The distribu- 
tions of source red-shift and total 2 — 10 keV counts (for the 
XIS-FI detector) are shown in Figure 1 and Figure 2, re- 
spectively. The AGN are predominantly local, with ~ 90% 
of the sample being located at a red-shift of z < 0.1, but also 
includes more distant objects such as PDS456 (z = 0.1840), 
1H 0419-577 (z = 0.1040), PKS 0558-504 (z = 0.1372) and 
PBCJ0839. 7-1214 (z = 0.198). The gravitationally lensed 
quasar APM 08279+5545 is by far the most distant object 
in the sample with a red-shift of z = 3.91. As shown in Fig- 
ure 2 the total 2 — lOkeV counts for the fitted spectra span 
two orders of magnitude but are approximately normally 
distributed about the peak and mean number of 10 5 counts. 
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Table 1. Source classifications 



Source classification 


Number 


Sy 1.0 - 1.2 


17 


Sy 1.5 


11 


Sy 1.8 - 1.9 


6 


NLSyl 


6 


BLRG 


6 


QSO 


5 


Total 


51 



3 DATA REDUCTION 
3.1 XIS Reduction 

There are four X-ray Imaging Spectrometers (XIS) CCD 
cameras aboard Suzaku which cover the energy range of 0.3- 
12.0 keV. The XISO, XIS 2 and XIS 3 are front-illuminated, 
while the XIS 1 is back-illuminated and provides superior 
sensitivity below ~ 1 keV but a lower effective area and a 
higher background at harder X-ray energies. The XIS 2 suf- 
fered a charge leak in November 2006 3 and is therefore only 
available in observations taken before this date. All spectra 
were reduced using v6.11 of the HE AS oft 4 suite following 
the process outlined in the Suzaku ABC Data Reduction 
Guide 5 . Events were selected using ASCA X-ray grades 0, 
2, 3, 4, and 6, adopting the standard screening criteria such 
that data were excluded if taken: (1) within 436 s of passage 
through the South Atlantic Anomaly (SAA), (2) within an 
Earth elevation angle (ELV) < 5°, and/or (3) with Earth 
day-time elevation angles < 20°. Hot and flickering pix- 
els were removed from the XIS images using the CLEANSIS 
script. Source spectra were extracted from within circular 
regions of radius 1.5' < r < 3.0' to ensure good coverage 
of the source events, while background spectra were typi- 
cally extracted from offset circles of the same radius with 
care taken to avoid the chip corners containing the Fe 55 
calibration sources. For intrinsically faint sources the back- 
ground spectra were extracted from circles larger than those 
of the source and the ratio of source/background area was 
accounted for with an appropriate background scaling fac- 
tor. Redistribution matrices and ancillary response files for 
each observation were generated using the tasks xisrmfgen 
and xissimarfgen, respectively. 

Where possible spectra obtained from the front- 
illuminated XIS detectors where combined into a single 
source spectrum (hereafter XIS-FI) using mathpha in or- 
der to maximize signal-to-noise (S/N) in the FeK band. 
Most of the observations in the sample have data for at least 
the XIS and XIS 3 with a further 22 also having data for 
the XIS 2 (see Table Al). For SWIFT J2127.4+5654 (OB- 
SID 702122010) only data for the front-illuminated XIS 3 
camera is available because the XIS was not operational 
during the observation. Spectra for the back-illuminated 
XIS 1 (hereafter XIS-BI) were reduced using the same pro- 
cess as outlined above and were analysed simultaneously as 
a separate input spectrum. All XIS spectra were grouped to 

3 http://heasarc.gsfc.nasa.gov/docs/suzaku/news/xis2.html 

4 http://heasarc.nasa.gov/lhcasoft/ 

5 ftp://legacy.gsfc.nasa.gov/suzaku/doc/suzaku_abc_guidc.pdf 



have a minimum of 50 counts per energy bin to enable the 
use of the \ 2 fit statistic. Net XIS exposures and the total 
background subtracted 2—10 keV count rates (in the source 
rest-frame), for both the co-added XIS-FI and the XIS-BI 
detectors, are listed in Table Al. 

3.2 HXD/PIN Reduction 

HXD/PIN spectra were also reduced using the method out- 
lined in the Suzaku Data Reduction guide and, again, pro- 
cessed according to the screening criteria described previ- 
ously. As the HXD/PIN is a collimating rather than an 
imaging instrument the contribution of both the instrumen- 
tal Non X-ray Background (NXB) and the Cosmic X-ray 
Background (CXB) need to be independently accounted for 
when estimating the total background. The NXB was gen- 
erated using the appropriate response and tuned event files 
(background model D) for each observation. Common good 
time intervals (GTIs) in the event and background files were 
selected using mgtime and spectral files were extracted us- 
ing xselect. PIN detector dead-time was accounted for us- 
ing the hxddtcor task and the NXB background exposure 
was increased by x 10 to reduce the effects of photon noise. 
The CXB contribution was simulated using the form of 
Boldt (1987), combined with the NXB to form a total back- 
ground spectrum using mathpha, which was subsequently 
subtracted from the source spectrum within xspec. 

All background subtracted PIN spectra were binned to 
at least the 3cr level above background (typically > 5a) to 
enable the use of x 2 fit statistic. Hard X-ray faint and/or 
high background observations where source count rates were 
< 4% of the PIN total were not considered in our analysis. 
Only APM 08279+5255 and PDS456, amounting to a total 
of five observations, meet this criteria (see Table Al). As 
with the XIS data, the final PIN exposures and the total 
background subtracted rest-frame 15 — 50keV counts are 
listed in Table Al. 



4 ANALYSIS 

4.1 Spectral Fitting 

A detailed broadband spectral analysis of all sources was 
conducted using version 12.6.0q of the xspec spectral fit- 
ting package. All spectral models are modified by the ap- 
propriate column density of Galactic absorption using val- 
ues taken from the Leiden/ Argentine/Bonn Survey of Galac- 
tic HI (Kalberla et al. 2005) which were obtained from the 
on-line version of the Nh FTOOL 6 . The values of Galactic 
absorption for each source are listed in Table Al. XIS-FI 
data were typically considered between 0.6 — 10.0 keV, while 
XIS-BI were only included between 0.6 — 5.0 keV due to de- 
creasing S/N in the FeK band which could hamper absorp- 
tion line detection. All XIS spectra were excluded between 
1.6 — 2.1 keV due to uncertainties in calibration associated 
with the SiK detector edge. Where available, PIN data were 
included to cover at least the 15 — 40keV energy range. 

A constant multiplicative factor was included in all 
models to account for the XIS/PIN cross-normalization 

6 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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whose value depends not only on the nominal pointing of the 
observation, but also the version of the Suzaku pipeline with 
which the data has been processed. The current XIS/PIN 
ratios suitable for version 2 processed data are 1.16(1.18) 
for the XIS(HXD) nominal pointing positions 7 . However, 
the cross-normalization was up to ~ 5 — 6% lower for data 
processed with version 1 of the pipeline, corresponding to 
a XIS/PIN ratio of 1.09(1. 13) 7 in the XIS(PIN) nominal 
pointing positions. While the difference is only small the 
additional uncertainty can have a considerable effect on the 
continuum parameters at E > lOkeV, particularly in hard 
X-ray bright sources with a high PIN count rate where the 
model can become driven by the hard X-ray band. To ac- 
count for any systematic effects associated with the instru- 
mental cross-calibration we therefore allow the constant pa- 
rameter in each model to vary ±5% about the values sug- 
gested by the Suzaku team to take into account any system- 
atic errors. 

There are 20 AGN in the sample which have 
two or more Suzaku observations. In APM 08279+5255, 
IC4329A, MCG -6-30-15, Mrk841, NGC5506, NGC 5548 
and PKS 0558-504 the different observations are similar in 
spectral shape which allows them to be co-added using the 
appropriate relative weighting factors to take into account 
differences in individual exposures; with the resultant time- 
averaged spectra for these sources being used in all subse- 
quent analyses. 1H 0419-577 and NGC 2992 both had ob- 
servations which were taken in different Suzaku nominal 
pointing positions which could influence co-adding. How- 
ever, effective area was always consistent to within ±10% 
and so the spectra were still co-added using the mean of 
the response files. Any additional systematic uncertainty 
introduced was adequately accounted for by the variable 
XIS/PIN cross-normalisation. In 3C 120 and Mrk509, which 
both showed notable variability between observations, spec- 
tra were jointly fit depending on the extent of their spectral 
variability between epochs. In 3C 120 we followed the analy- 
sis method of Kataoka et al. (2007) and Tombesi et al. (2010) 
by co-adding OBSIDs 700001020, 70001030 and 70001040, 
which are all of a similar spectral shape and flux level, and 
jointly fitted them with OBSID 700001010 which has a more 
prominent underlying soft-excess. In Mrk509, which is well 
known for having a strong soft-excess (Mehdipour et al. 
2011), OBSID 701093010, the stacked OBSIDs 701093020, 
701093030, 701093040, and OBSID 705025010 were fitted 
simultaneously to account for the observed variability in the 
soft X-ray band. 

The remaining nine sources with more than one ob- 
servation (3C111, Fairall9, Mrk766, NGC 1365, NGC 3227, 
NGC 3516, NGC 3783, NGC 4051 and PDS456) showed 
strong spectral variability in both the shape of the spectrum 
and/or drastic changes in flux state which made co-adding 
impractical. In these sources the available spectra were fit- 
ted simultaneously and a model was constructed to describe 
the observed spectral variability with as few additional free 
parameters as possible. In NGC 1365, OBSIDs 702047010 
and 705031010 are dominated by very deep FeK absorp- 
tion lines. These lines are not present in OBSID 705031020, 

7 http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzak' 
2007-ll.pdf 



possibly due to the source dropping into a quasi-scattering- 
dominated state (see Maiolino et al. 2010 for details of the 
variability patterns in this source). For simplicity we there- 
fore only simultaneously fitted OBSID 702047010 and OB- 
SID 705031010 in NGC 1365, and fit OBSID 705031020 sep- 
arately. 

The x 2 minimization technique is used throughout this 
work; with all statistical errors quoted to the 90% confidence 
level (Ax 2 = 2.71 for one parameter of interest). Where the 
statistical significance of components is quoted in terms of 
a Ax 2 value the component in question has been removed 
from the model and the data refitted to ensure that the order 
in which components are added to the model has no influ- 
ence on the derived statistics. When referring to statistical 
changes to a fit a positive Ax denotes a worsening of the 
fit, while a negative Ax 2 indicates a statistical improvement. 
Positive outflow velocities correspond to a net blue-shift rel- 
ative to the systematic of the host galaxy, while a negative 
velocity indicates a net red-shift. 

4.2 Model Construction 

All spectra were first fitted with a power-law modified solely 
by Galactic absorption. Additional components were added 
to the model and retained provided their significance ex- 
ceeded the > 99% confidence level by the F-test. Any emis- 
sion lines in the soft X-ray band were fitted with narrow 
Gaussians (er ~ 5eV). In sources with multiple spectra we 
initially found a broadband fit to the observation with the 
highest flux; later observations were then added sequentially 
and the continuum/absorption parameters were allowed to 
to vary independently until a simultaneous fit to all spectra 
had been achieved using as few free parameters as possi- 
ble. Once a statistically acceptable fit to the broad-band 
spectrum of each source had been found, we fitted any nec- 
essary emission components and then performed a system- 
atic and methodical search for Fe K absorption lines between 
5 — lOkeV using energy-intensity plane contour-plots (Sec- 
tion 4.3) and detailed Monte Carlo simulations (Section 4.4). 

The modelling approach for individual spectral compo- 
nents are outlined below, with the continuum parameters 
for each source being noted in Table Dl and those for the 
warm absorber being listed in Tables D2 and D3 for single- 
and multi-epoch spectra, respectively. A description of the 
underlying modelling assumptions and associated caveats is 
presented in Appendix C. 

4-2.1 Warm absorption I: Fully covering 

Depending on the properties of the intervening material 
(such as its ionisation state and column density) absorp- 
tion by circumnuclear material can add significant spectral 
curvature to the observed X-ray spectrum and can therefore 
have a direct effect on the continuum and line parameters 
measured in broadband models. In this work we model warm 
absorption components using a suite of xstar (v. 2.2. lbc) ta- 
bles which are all generated assuming input values which are 
'typical' for local Seyfert galaxies 8 (e.g., Til). The resultant 

® Absorption grids arc described by: an illuminating photon in- 
dex of T = 2.1, a gas density of n = 10 10 cm~ 3 , a micro-turbulent 
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xstar tables cover a wide range of parameter space in terms 
of column densities [10 18 < log(W H /cm -2 ) < 10 24 ], and ion- 
isation parameter [—3 < log(£/ergcms~ ) < 6] which makes 
them well suited for accounting for all manner of warm ab- 
sorption. 

Fully-covering warm absorption zones are included in 
models where necessary to achieve a good fit to the soft X- 
ray band; in some cases more than one absorption zone is 
needed. In these cases the column density and ionisation pa- 
rameter of each zone is allowed to vary independently, and 
represents an absorption geometry which consists of multi- 
ple layers of gas. At the energy resolution of the XIS CCDs 
the bound-bound absorption lines required to constrain the 
outflow velocities of individual soft X-ray absorption com- 
ponents are unresolved and all absorption zones are there- 
fore fitted as stationary in the source rest-frame (i.e., fixed 
outflow velocities of v ou t = Okms -1 ). Allowing the outflow 
velocity of the soft X-ray absorber to vary always has a negli- 
gible effect on the reported Fe K absorption line parameters. 

4-2.2 Warm absorption II: Partial covering 

We also consider the possibility that the sight-line to a 
source is partially covered. In this absorption geometry a 
fraction / cov < 1 of the source flux is absorbed with the re- 
maining 1 — / cov leaking through the absorption layer. For 
simplicity we account for partially-covering absorption lay- 
ers using a customized version of the zxipcf model 9 which 
models the partially-covering absorption by partially-ionised 
gas (see Reeves et al. 2008a) without needing to use compli- 
cated nested power-laws and xstar tables. As in zxipcf the 
free parameters in our customized table model are column 
density (A/h), ionisation parameter (logf), covering fraction 
(/cov) and red-shift relative to the observer (z). The model 
is based on the same tables as discussed in Section 4.2.1 and 
therefore has a slightly lower turbulent velocity than zxipcf, 
at Utmb = 100 km s -1 , but covers the same parameter space 
in terms of column density and ionisation parameter. 

Partially-covering absorption can have a strong effect 
on the observed continuum with moderate column densities 
of material (Ah ~ 10 23 cm~ 2 ) adding considerable spec- 
tral curvature at E < lOkeV (Reeves et al. 2004; Risaliti 
et al. 2005; Braito et al. 2011; Turner et al. 2011). Bound- 
free transitions in similar column density material can also 
fit broad residual emission profiles in the FeK band (In- 
oue & Matsumoto 2003; Miller et al. 2008, 2009; Tatum 
et al. 2012), and partial-covering by Compton-thick mate- 
rial (A^h > 10 24 cm~ 2 ) can reduce the observed emission be- 
low 10 keV with the true intrinsic continuum only becoming 
apparent at higher energies as a 'hard excess' of emission 
relative to that expected from standard reflection models 
(Reeves et al. 2009; Turner et al. 2009; Risaliti et al. 2009; 
Tatum et al. 2012). 

We include partial-covering absorbers in our models if 
and when they are required by the data at the Pp > 99% 

velocity of i>turb = 100 km s -1 , and an integrated model luminos- 
ity of L = lO^ergs- 1 between 1 - 1000 Rydbergs. 
9 The standard zxipcf uses a specific grid with the following 
parameters: T = 2.2, n = 10 10 cm~ 3 , v tU ib = 200kms _1 , L = 
10 44 ergs" 1 



confidence level and a satisfactory fit to the data could not 
be achieved using solely fully-covering absorption; several 
sources appear to need multiple partially-covering absorp- 
tion zones which suggests the presence of a clumpy stratified 
absorber along the line of sight. Again, all absorber param- 
eters are listed in Tables D2 and D3 in Appendix D. We 
note that in some circumstances the need of a high column 
density partially covering component can be contingent on 
the means with which the underlying reflection component 
is modelled leading to some model degeneracies. However, 
as we discuss in Appendix C, regardless of whether the hard 
X-ray data is fitted with reflection or not, partial-covering 
has little measurable effect on the parameters measured for 
any highly-ionised absorption line systems. 

4.2.3 The Soft- excess 

Relative to the low energy extrapolation of the power-law 
continuum in the 2 — lOkeV band the X-ray spectrum of 
AGN often show a smooth increase in emitted flux below 
~ IkeV (Turner & Pounds 1988; Porquet et al. 2004). The 
thermal temperature of this 'soft-excess' suggests it is un- 
likely to be the direct emission from a standard accretion 
disc without additional reprocessing (Sobolewska & Done 
2007; Done et al. 2012). Alternative explanations posit that 
the soft-excess may be due to an increase in optical depth 
associated with circumnuclear O VII- VIII and Fe L-shell tran- 
sitions at E < 0.7keV which can enhance either the trans- 
mitted or reflected flux along the sight-line through smeared 
absorption (e.g., Chevallier et al. 2006; Sobolewska & Done 
2007; Done &: Nayakshin 2007) or blurred reflection (e.g., 
Crummy et al. 2006; Nardini et al. 2011; Brenneman et al. 
2011; Nardini et al. 2012) effects. Furthermore, in some 
sources (e.g., Mrk766, Miller et al. 2007; MCG -6-30-15, 
Miller et al. 2008) the 'excess' could simply be a product of 
complex absorption and just be the manifestation of opacity 
around ~ 1 — 2 keV. 

Regardless of the true physical origin we take a purely 
phenomenological approach when fitting the soft-excess in 
this work, and predominantly use the bbody model which 
represents the emission from a constant temperature black- 
body. While not necessarily physically motivated, modelling 
the soft-excess in this manner offers a simple parametrisa- 
tion of any 'excess' soft X-ray emission which is sufficient to 
get a good handle on the underlying continuum parameters. 
For completeness we investigate the effects that other ways 
of modelling the soft-excess can have on any Fe K-band ab- 
sorption lines in Appendix C. Roughly half of the sources in 
the sample (24/51; ~ 47%) show evidence for a soft-excess, 
of which 22/24 are fit with a bbody component. In 3C 120 
and Mrk509 the soft-excess is very broad and extends be- 
yond that which can be fitted with a simple blackbody. In 
these sources we instead fit the excess with a second power- 
law with a softer photon-index. 

4-2.4 Lowly ionised reflection 

Cold reflection from large column densities of neutral or 
lowly ionised material outside of the sight-line can have 
a strong influence on the observed X-ray spectrum. The 
strongest observational characteristics of such reflection in- 
clude the Compton-reflection hump at ~ 30 — 40 keV and 
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the almost ubiquitous Fe Ka and Fe K/3 fluorescence lines at 
~ 6.4 keV and ~ 7.06 keV (Nandra et al. 2007; Shu et al. 
2010), respectively. Compton down-scattering of both Ka 
line photons and high energy continuum emission also gives 
rise to a 'Compton shoulder' at ~ 6.2 keV (e.g., Matt 2002; 
Yaqoob & Murphy 2011) and resonant line emission in the 
soft X-ray band (e.g., Ross & Fabian 2005; Garcia & Kali- 
man 2010) which can further complicate the emergent spec- 
trum. 

Naturally, owing to the important effect it can have 
on the observed X-ray spectrum there are numerous mod- 
els available for modelling the reflection component (e.g., 
pexrav/pexriv, Magdziarz & Zdziarski 1995; pexmon, Nan- 
dra et al. 2007; reflionx, Ross & Fabian 2005; xillver 10 , 
Garcia & Kallman 2010; MYTorus 11 , Murphy & Yaqoob 
2009). In this work we use a combination of reflionx and 
pexrav, both of which are publicly available and extensively 
used in the literature. The primary reason for this it that 
because reflionx interpolates the observed reflection spec- 
trum from a pre-generated grid of table models it is signif- 
icantly faster at fitting spectra than, for example, pexrav, 
pexriv or pexmon, which analytically calculate the Comp- 
ton reflection spectrum on the fly at each step of the fitting 
process, or MYTorus which requires the model to be tailored 
for each individual source. Using reflionx as our primary 
means of fitting the reflection spectrum ensures the least 
time-consuming Monte Carlo simulations which is impor- 
tant when dealing with a large sample of objects such as 
that considered in this work. Secondly, when fitting reflec- 
tion continuum without the simultaneous constraint of the 
FeKa line the reflection fraction, R, reported by pexrav 
(and pexriv) can become degenerate with the photon-index 
of the primary power-law, with a hardening reflection com- 
ponent compensated for by a softer primary continuum. By 
simultaneously fitting the reflection continuum, the FeKa 
line and the soft X-ray resonance lines, reflionx is able to 
overcome these modelling degeneracies which leads to a con- 
fident constraint on the contribution of the reflection con- 
tinuum to the observed spectrum, reflionx also has the 
additional advantage of allowing the ionisation state of the 
reflector, £, to be a free parameter which enables it to model 
changes in the FeKa emitted flux and FeK-shell edge pro- 
file associated with the reflectors ionisation state. We stress, 
however, that equivalent results are found for the detected 
FeK absorbers if pexmon is used instead, but with the resul- 
tant Monte Carlo simulations taking significantly longer to 
complete which effectively prohibits its uniform use through- 
out the sample. 

We initially fitted all sources with pexrav to de- 
termine the parameters of the FeKa and FeK/3 lines, 
pexrav was then replaced with reflionx, and a system- 
atic search for additional atomic features in the Fe K band 
was conducted. A total of 11 sources have best-fitting re- 
flector Fe abundances which are non-solar, with 4 requir- 
ing a slight over-abundance (MCG +8-11-11, NCG4593, 
NCG7213, NCG7469) and 7 with an under-abundance 

10 The xillver reflection model is not currently available in the 
public domain and is therefore not considered for use in this work. 

11 The MYTorus model and documentation are publicly available 
at: www.mytorus.com 



(4C +74.26, Fairall9, IC4329A, IGR J16185-5928, Mrk335, 
Mrk359, Mrk841). These abundances are most likely a by- 
product of reflionx assuming a face-on reflection geometry, 
and is a caveat which is discussed in greater detail in see Ap- 
pendix C. PDS 456 and APM 08279+5255 are not fitted with 
a reflection component because they lack sufficient counts in 
the HXD/PIN.. 

4.3 Searching for Fe K absorption 

Once a statistically acceptable fit to the broadband contin- 
uum had been found, i.e., including all necessary absorption 
regions, soft X-ray emission lines and continuum compo- 
nents, we performed a thorough search for additional spec- 
tral features between 5 — 10 keV. The method consists of 
inspecting the \A\ | deviations from the best-fit continuum 
model using inverted contour plots of the energy-intensity 
plane in the Fe K band. The method of calculating the con- 
tour plots was adapted from the method outlined in T10A, 
and was carried out as follows: 

(i) an unresolved (a = 10 eV) Gaussian was stepped 
across the entire 5—10 keV energy band of the base- 
line continuum model in 25 eV intervals, with normal- 
ization allowed to adopt both positive and negative 
values to probe for spectral lines in both emission and 
absorption. All of the other spectral components were 
allowed to be free; 

(ii) after each step the A\ 2 deviation was recorded gener- 
ating a x 2 distribution of the entire Fe K band relative 
to the baseline continuum model; 

(iii) confidence contours for the grid of \ 2 values were 
plotted according to A\ 2 deviations of —2.3, —4.61, 
—9.21, —13.82, —18.42 from the baseline model, which 
correspond to confidence intervals of 68%, 90%, 99%, 
99.9% and 99.99%, respectively. A confidence contour 
corresponding to a A\ 2 = +0.5 worse fit is also plot- 
ted which is intended to denote an approximate level 
for the continuum. 

The energy-intensity contour plots produced with this 
method provide a powerful means of searching for additional 
emission or absorption components present in the Fe K band 
while also visually assessing their energy, intensity and rough 
statistical requirement relative to the underlying continuum 
model. 

There are a number of atomic features between 6.0 — 
8.0 keV which can complicate the detection of absorption 
line systems. Such features include ionised emission lines 
from Fexxv and Fexxvi expected at ~ 6.63 — 6.7keV and 
~ 6.97 keV, respectively, the FeK-shell edge complex at 
~ 7.1 keV, and any broad Fe line profile due to GR or scat- 
tering effects. It is important to account for these emission 
residua prior to searching for absorption lines, particularly 
in sources which show evidence for a broad residual as the 
broadness of the profile can effectively mask the presence of 
low velocity absorption systems in the raw data. For a given 
continuum model we searched for highly-ionised absorption 
lines in the Fe K-band using the following steps: 

(i) we first generated an energy-intensity contour plot us- 
ing the method outlined above; 

(ii) we then inspected the contour plot to determine 
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Ark 120 (OBSID: 702014010) 




7 8 9 

Rest-frame Energy (keV) 

Figure 3. Top panel: Ratio plot showing the residuals remain- 
ing in the FeK band of Ark 120 (OBSID 702014010) once all 
atomic lines have been removed and reflection fitted with pexrav; 
Middle panel: confidence contours showing deviations from the 
best-fit model when the FeKo and FeK/3 lines have been fitted 
with reflionx and a narrow Gaussian respectively. The closed 
significance contours correspond (from outer to inner) to A% 2 
improvements of -2.3 (P F = 68%), -4.61 (P F = 90%), -9.21 
(P F = 99%), -13.82 (P F = 99.9%) and -18.42 (P F = 99.99%) 
relative to the best-fit continuum model are shown in red, green, 
blue, cyan and magenta, respectively. The magenta contours in 
the middle panel indicate that the broad emission residual is sig- 
nificant at the > 99.99% level and that an additional component 
is required in the model; Bottom panel: the remaining confidence 
contours once the broad profile has been fitted. The broad residua 
are no longer detected but an additional Fexxvi emission line is 
present at E ~ 6.97 keV which is significant at > 99%. In all pan- 
els the dashed vertical lines indicate the expected rest-frame ener- 
gies of (from left to right) the FeKct fluorescence line, Fexxv Hect 
and Fexxvi Lya, respectively. Colour version available online. 



whether there were any intense ionised emission and 
absorption lines present in the data with confidence 
contours of > 99%; 
(iii) where there was evidence for a broad emission residual 
at an F-test significance of at least P F > 99% they 
were fitted with either a broadened Gaussian (with a- 



width a free parameter) or a diskline depending on 
the asymmetry of the observed profile, and a second 
intermediate contour plot was generated to determine 
whether any further components were needed by the 
model. As before, all other model parameters were 
allowed to vary freely during this process. If there was 
no evidence for a broad profile we did not generate 
an additional contour plot and instead moved directly 
onto the next step; 

(iv) where narrow emission profiles were detected with a 
resolved confidence contour of > 99% they were fitted 
with unresolved (a = 10 eV fixed) Gaussian profiles 
provided they were required by the data at P F > 99% 
(corresponding to A\ 2 > 9.21 for two parameters of 
interest); 

(v) once all emission profiles had been accounted for we 
again checked for the presence of any blue-shifted Fe K 
absorption lines at E > 6.7 keV. If no absorption lines 
were present we ended the search at this step, reported 
the best-fit continuum and FeK emission line param- 
eters in the relevant tables of Appendix C, and moved 
onto the next observation in the sample; 

(vi) otherwise, where absorption troughs were clearly de- 
tected with a confidence contour of > 99% we 
parametrised the line(s) using inverted Gaussians with 
cr-width initially fixed at either 10 eV, 30 eV or 100 eV 
depending on which provided the greatest improve- 
ment to Ax 2 - Note that while initially fixed, the line 
widths were allowed to vary where appropriate lead- 
ing to four sources with resolved absorption profiles. 
The key parameters of all detected absorption lines 
are reported in Table 2. 

This process was carried out on each individually fitted spec- 
trum in the sample; including those which were included in 
simultaneous fits. In Figures 3 and 4 we show examples of 
this process applied to Ark 120 (OBSID 702014010), which 
shows evidence for a broad asymmetric emission profile and 
no absorption lines, and Mrk 766 (OBSID 701035020), which 
is dominated by Fexxv Hea and Fexxvi Lya absorption. 
The top panel of both figures show the residuals which re- 
main in the FeK band when all atomic lines have been re- 
moved from the best-fit continuum model and the reflec- 
tion component is fitted with pexrav to highlight the pres- 
ence of the neutral Fe Ka/Fe K/3 fluorescence lines. The con- 
tour plots show the significances of the remaining residu- 
als when the FeKa and FeK/3 lines have been fitted with 
reflionx and a narrow Gaussian line. The continuous outer 
contour corresponds to the A\ 2 = +0.5 residual as men- 
tioned previously. From outer to inner the closed signifi- 
cance contours corresponding to Ax 2 improvements of —2.3 
(P F = 68%), -4.61 (P F = 90%), -9.21 (P F = 99%), -13.82 
(P F = 99.9%) and -18.42 (P F = 99.99%) relative to the 
best-fit continuum model are shown in red, green, blue, cyan 
and magenta, respectively. 

12 For simplicity, the diskline was fitted with an assumed rest- 
frame energy of E = 6.4 keV and a typical emissivity profile of 
q = —2.5 (e.g., Patrick et al. 2011). The inner emission radius 
(Pin), outer radius (Pout) and disk inclination (i), were left free 
to vary. 
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Mrk766 (OBSID: 701035020) 
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Figure 4. Top Panel: As in the case of Figure 3, but this time 
for Mrk766 (OBSID 701035020); Bottom Panel: As in the corre- 
sponding panel of Figure 3. There are no ionised emission lines 
required in Mrk 766 however two highly significant (each required 
at Pp > 99.99%) absorption profiles are clearly detected. The en- 
ergy of these lines is consistent with Fexxv Hect and Fexxvi Lyct, 
respectively, outflowing at v ou t ~ 6000 km s — 1 ). Colour version 
available online. 



There are clear residual profiles present in both sources; 
particularly Ark 120 where the positive asymmetric profile 
and additional Fexxvi emission is apparent in the middle 
and bottom panels of Figure 3, respectively. The two ab- 
sorption profiles in Mrk 766 - which have previously been 
identified as Fexxv Hea and Fexxvi Lycv by Miller et al. 
(2007) - are statistically distinguishable at > 99.99% con- 
fidence and are self-consistently fitted with a single highly- 
ionised region of photoionised absorption (see Section 4.6). 
Energy-intensity contour plots for all fitted spectra in this 
sample are included in Figures Bl and B2 in Appendix B. 
We note that the A^ 2 of individual lines and their corre- 
sponding significances according to the F-Test as listed in 
Table 2 are taken from a direct spectral fitting of the line 
profiles themselves, and are not determined directly from 
the contour plots. 

It is also important to note that while we use the F-test 
as a rough initial gauger of the significance of any line-like 
profiles in the FeK band we do not claim to robustly de- 
tect any absorption lines based solely on their measured A\ 2 
and corresponding F-Test significances. Indeed, as has been 
pointed out by numerous authors in the literature (e.g., Pro- 
tassov et al. 2002; Porquet et al. 2004; Markowitz et al. 2006; 
T10A) the F-test alone might not be an adequate statistical 
test to determine the detection significance of atomic lines 
in complex spectral models as it neither takes into account 
the number of energy resolution bins over a given energy 



range, or the expected energy of a given atomic line. In cases 
where there is no a priori justification for expecting a spec- 
tral line at a particular energy and the line search is done 
over what is essentially an arbitrary energy range, as is often 
the case when dealing with strongly blue-shifted absorption 
lines, the F-test can somewhat over-predict the detection 
probability and when compared to extensive simulations. 
For these reasons all suspected absorption lines which have 
an F-test significance of Pf > 99% were followed up with 
extensive Monte Carlo simulations which allows their detec- 
tion significance to be assessed against random fluctuations 
and photon noise in the spectrum. 

4.4 Monte Carlo simulations 

Such Monte Carlo simulations have been used extensively 
in the literature to overcome the limitations of the F-test 
(e.g., Porquet et al. 2004; Markowitz et al. 2006; Miniutti & 
Fabian 2006; T10A; T10B) and enable the statistical signif- 
icance of a spectral line to be robustly determined indepen- 
dently of spectral noise and detector effects. The method 
of Monte Carlo simulation we used follows the same pro- 
cess which was first outlined by Porquet et al. (2004), and 
is almost identical to that used by T10A. The process was 
carried out with the following steps: 

(i) from the null hypothesis model (i.e., the broadband 
continuum model with all narrow absorption lines 
in the FeK band removed) we simulated both XIS- 
FI and XIS-BI spectra using the fakeit command in 
XSPEC and subtracted the appropriate background 
files. The simulated spectra had the same exposure as 
the original data and used the same spectral response 
files; 

(ii) we then fitted the simulated XIS-FI(XIS-BI) spectra 
between 0.6 - 10.0keV(0.6 - 5.0 keV) with the null 
hypothesis model to produce a new and refined null 
hypothesis which takes into account the uncertainty 
in the null hypothesis model itself. All continuum pa- 
rameters bar the photon-index of the primary power- 
law and its normalization, and the normalization of 
any bbody component, were frozen to their best-fit pa- 
rameters taken from the real data to prevent degen- 
eracies between model components during re-fitting. 
Any broad profile at FeK had its width fixed to the 
best-fit value but was allowed to vary in both centroid 
energy and normalization; 

(iii) from the refined null hypothesis model we generated 
a second set of simulated XIS-FI and XIS-BI spec- 
tra and subtracted the appropriate background files. 
These second simulated spectra were then fitted with 
the null hypothesis model and the resultant Xnuii value 
was recorded; 

(iv) an unresolved (<r = 10 eV) Gaussian was added to the 
model at 5 keV in the source rest-frame with intensity 
initially set to zero but left free to vary between both 
positive and negative values to probe for both emis- 
sion and absorption lines. The Gaussian line was then 
sequentially stepped between 5 — 9.5 keV (rest-frame) 
in 25 eV increments. After each step the |Axnoisel was 
recorded relative to Xnuiii 

(v) this process was carried out for T = 1000 simu- 
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Table 2. Gaussian parameters for detected FeK band absorption lines. Notes: (1) Source name; (2) Suzaku observation ID; (3) 
Absorption line identification; (4) Measured line energy in the source rest-frame, in units of kcV; (5) Measured line energy-width, in 
units of eV. Unresolved lines were fit with widths fixed at 10 eV, 30 eV or 100 eV depending upon which yielded the best statistical 
improvement to the fit; (6) Line intensity, in units of X 10 — 6 photons cm -2 s — 1 ; (7) Line equivalent width, in units of eV; (8) Change 
in fit statistic (and degrees of freedom) when line is removed from the best-fit model; (9) Line significance according to the F-test, 
in per cent; (10) Monte Carlo significance of individual lines, in per cent; (11) Monte Carlo significance of observing a pair of lines 
with energy separation corresponding to a Fexxv-xxvi pair (see text for details), in per cent. Absorption complexes detected at 
Pmc ^ 99% confidence are listed in bold. 
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lated spectra per observation yielding a distribution of 
|Axnoisel under the null hypothesis which essentially 
maps the statistical significance of any deviations from 
the null hypothesis model which are solely due to ran- 
dom photon noise in the spectrum; 
(vi) the measured significance of the line in the real spec- 
trum |Axnnel was then compared to the values in the 
|AXnoisel distribution to assess how many simulated 



spectra had a random fluctuation with a detection 
significance over this threshold value. If N simulated 
spectra have |Ax no isel ^ I^Xnnel then the estimated 
confidence interval for the observed line from Monte 
Carlo simulations is then Pmc" = 1 — y ^ y nc j . More- 
over, if there are two absorption profiles consistent 
with a Fe xxv-xxvi pair we can infer the null proba- 
bility of both lines simultaneously being a false detec- 



© ? RAS, MNRAS 000, 1-?? 



The Suzaku view of highly-ionised outflows in AGN. 



11 



tion by multiplying the probabilities of each individual 
line. 

Monte Carlo significances for all absorption lines detected in 
this work are listed in column 10 of Table 2. Absorption com- 
plexes with a total Pmc =S 99% are conservatively identified 
as robustly detected, while those with 95% ^ -Pmc < 99% 
are only listed as marginal detections. The overall detection 
rate and global significance of Fe K-band absorption lines is 
further discussed in Section 5.1. 

There are two possible caveats associated with the 
above Monte Carlo process which both warrant further dis- 
cussion. First and foremost, the Monte Carlo simulations im- 
plicitly assume that the null-hypothesis model is the correct 
representation of the continuum in a given source, and there- 
fore the Monte Carlo probabilities do not account for the 
possibility of continuum mis-modelling. Given this possibil- 
ity it is important to note that we have attempted to achieve 
a statistically acceptable representation of the broad-band 
spectrum, so that no obvious broad residuals are present. 
Care has been taken in the FeK band in particular such 
that any broad residuals are minimised prior to searching 
for absorption lines, such that the reduced-x 2 is ~ 1.0 in 
all cases. Moreover, and as shown in Appendix B, we note 
the vast majority of suspected absorption residuals are man- 
ifested through discrete narrow dips in the spectrum rela- 
tive to the best-fit model, whereas any systematic residuals 
would usually be broader than the instrument resolution. 

The second caveat is associated with spectral complex- 
ity in the Fe K band which can complicate both line identifi- 
cation and spectral interpretation. There can be significant 
atomic complexity between ~ 5 — 7keV, e.g., the narrow 
FeKa and FeK/3 fluorescence lines, broad underlying FeK 
lines, and ionised Fexxv-xxvi emission lines, and the de- 
tection of absorption in this regime can depend strongly on 
how these features are modelled. There are three sources 
(4C +74.26, MGC-6-30- 15, SWIFT J2127.4+5654; see Fig- 
ure B2) in the sample in which both a broad underlying 
profile and at least one absorption line have been detected. 
In both 4C+74.26 and SWIFT J2127.4+5654 the absorption 
lines are detected at a high energies (i.e., E > 8 keV in the 
source rest-frame) and the effect of the broad Fe line on the 
line detection is negligible. In MCG -6-30-15 the lines have 
also been confirmed with other X-ray observatories (Young 
et al. 2005; Miller et al. 2009) which suggests that the pres- 
ence of a broad Fe line does not introduce any significant 
model systematics. Thus modelling of a broad emission line 
does not appear to effect the detections of Fe K absorption 
lines in these cases. 



4.5 Consistency Checks 

To further test the robustness of the absorption lines de- 
tected in the co-added XIS-FI spectra we performed a se- 
ries of consistency checks with the individual XIS detectors. 
If a line is detected in two (or more, if the XIS 2 is also 
present) detectors the line is very likely to be a real fea- 
ture intrinsic to the source spectra rather than an artefact 
of background subtraction, systematic noise or an associ- 
ated detector effect. The XIS-BI is not suitable for use as a 
consistency check as it tends to have a much lower signal-to- 
noise ratio above around 5 keV owing to its lower effective 
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Figure 5. Consistency check comparing the absorption line 
equivalent width measured with the XIS versus the XIS 3 along 
with the 90% error bars. The diagonal dashed line represents the 
location where the EW is equal in both detectors. The circle 
(black), diamond (blue) and square (red) markers show the EW 
of the Fexxv Hea, Fexxvi Lya and blended Fexxv-xxvi lines, 
respectively. For sources with lines detected in more than one 
observation the mean EW for each line in each detector and as- 
sociated errors are plotted. Where available this consistency check 
was also carried out with the XIS 2 which yielded results which 
are entirely consistent with those displayed here. Colour version 
available online. 



area and higher background. Where absorption has been 
detected at the Pmc =S 95% level we fitted each individ- 
ual background-subtracted XIS 0, XIS 2 (where present) and 
XIS 3 spectrum with the best-fitting model to the co-added 
XIS-FI spectrum. We then fitted a Gaussian absorption line, 
with rest-frame energy and normalization left free to vary, 
at the energy where the absorption line is detected in the 
XIS-FI spectrum, refitted the data, and noted the resul- 
tant line parameters for each XIS detector. We note that in 
SWIFT J2127.4+5654 this consistency check could not be 
conducted as only the XIS 3 spectrum was available. 

In all sources the residuals detected in the individual 
XIS spectra have Gaussian parameters consistent with those 
found for the co-added FI spectrum. In Figure 5 we show a 
comparison of absorption line equivalent width (EW) , along 
with the 90% error bars, as measured independently with the 
XIS and XIS 3 detectors. The dashed diagonal line repre- 
sents the position where the EW is equal in each of the 
detectors; where the EW is consistent at the 90% level the 
error bars will overlap this line. In this plot Fexxv Hea, 
Fe xxvi Lya and blended Fe xxv-xxvi lines are plotted sep- 
arately, even in the cases where the lines are found to be part 
of a Fe xxv pair. We find that the mean EWs measured with 
the XIS and XIS 3 detectors are consistent at the 90% level 
(i.e., the error bars cross the diagonal line in Figure 5) in all 
sources. Where possible this analysis was also carried out 
using the XIS 2 spectra versus both the XIS and XIS 3. 
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Again, in all cases the mean parameters of the lines were 
always consistent at the 90% level. The consistency check 
strongly suggests that the observed lines are real features 
and are not due to detector effects or due to background 
subtraction. 

4.6 Photoionisation modelling 

All absorption complexes detected at a Pmc 95% confi- 
dence level were then fitted with the xstar photoionisation 
code to probe the likely physical properties of the absorbing 
material, such as the column density Nn, ionisation param- 
eter £, and the red-shift of the absorber relative to the ob- 
server z a . Probing these parameters is important as not only 
do they allow the mean properties of the absorbing material 
to be determined but, through the use of simple geometric 
assumptions, they also permit an order-of-magnitude assess- 
ment of an absorbers likely radial distance from the ionis- 
ing source, the mass outflow rate, and their global energetic 
output. Detailed discussion regarding the the absorber kine- 
matics will be presented in a companion paper (Gofford et 
al., in prep). 

In this work the absorption lines were fitted with 
an xstar table that had an assumed illuminating contin- 
uum of r = 2.1 and a micro-turbulent velocity (« tur b) of 
lOOOkms -1 , which roughly corresponds to the full-width 
at half-maximum (FWHM) velocity width of a a = 10 eV 
Fexxvi Lya absorption line at a rest-frame energy of 
6.97 keV. In cases where the lines were resolved (i.e., a > 
10 eV) we used xstar tables with a «t ur b values which closer 
matched the measured FWHM velocity width of the ob- 
served profile we used; we found that tables with v tU ib equal 
to 1000 kms"\ 3000 kms" 1 , 5000 kms" 1 and 10000 kms" 1 
were sufficient to fit all systems in the sample. In 17/20 
sources with significantly detected absorption lines all of the 
parameters of the xstar table were allowed to vary freely. 
In APM 08279+5255 and the joint spectra fits to NGC 1365 
and PDS 456 some parameters were tied to prevent degen- 
eracies between the column density and ionisation parameter 
(see Table 3). The absorber red-shift (z ) as measured from 
the spectrum using xstar is given in the observer frame. 
This is related to the intrinsic absorber red-shift in the 
source rest-frame (z a ) and the cosmological red-shift of the 
source (z c ) through the relation: (1 + z a ) = (1 + z a )(l + z c ). 
From this, the intrinsic velocity of the outflow relative to 
the source (iw) can then be calculated from the relativistic 
Doppler formula which ensures that the relativistic effects 
associated with both high red-shift sources and high- velocity 
outflows are correctly taken into account when inferring ab- 
sorber outflow velocities relative to the source rest-frame. 

4-6.1 Line identifications 

Before discussing the results of the xstar fitting it is im- 
portant that the a priori assumption that the absorption 
lines detected at E > 6.6 keV are due to the velocity 
shifted resonance lines of Fexxv and Fexxvi is justified. 
Indeed, while the Ka transitions of Fexxv and Fexxvi, 
which are expected at mean rest-frame energies of 6.697 keV 
and 6.966 keV, respectively, are expected to be the strongest 
lines in the « 6.5 — 7.0 keV energy interval, there are several 



other atomic features at higher energies which may compli- 
cate the identification of blue-shifted absorption systems at 
E > 7kcV. For example, the K-shell edges from the various 
ionised species of Fe are found above 7keV, with energies 
ranging from 7.1 keV for neutral Fe up to 9.3 keV for Fe xxvi. 
In the case of low-moderate ionisation Fe the K-edge is ac- 
companied by higher order resonance line structure which 
can give the edge a subtle curved profile in CCD spectra 
rather than it simply being an abrupt drop in flux (Kali- 
man et al. 2004). Furthermore, given its proximity to the 
rest-frame energy of the Fe xxv-xxvi transitions the neutral 
edge at 7.1 keV could have an influence on the detection of 
lower velocity systems. 

Of the sources included in the sample CBS 126, 
Mrk766 (OBSIDs 701035010, 701035020), NGC 3227 (OB- 
SIDs 703022010, 703022020, 703022030), NGC 4051 (OB- 
SIDs 700004010, 703023010) and NGC 4151 all have at least 
one significantly detected absorption line at a rest-frame en- 
ergy which is consistent at the 90% level with the neutral 
Fe K-shell edge from the reflection component. The absorp- 
tion in Mrk766, NGC 3227, and NGC 4051 is manifested by 
two lines with a common velocity shift equal to that ex- 
pected for Fexxv Hea and Fexxvi Lya and are therefore 
unlikely to be affected by the presence of an edge. However, 
the absorption in CBS 126 and NGC 4151, which only com- 
prises a single detected profile, could possibly be affected. 
Even so, given that the K-shell edge structure is already 
self-consistently accounted for by the reflionx model and 
in xstar, suggests that the residual profiles detected near the 
Fe K-shell edge in these sources is a real additional compo- 
nent rather than the residuals left by an inadequately fitted 
edge. Alternatively, the residuals could be due to a partially 
covering absorber with a low outflow velocity. We investigate 
this possibility further in Appendix C. Finally, we note that 
while NGC 5506 shows evidence for both a highly ionised 
Fexxv emission line at ~ 6.63 keV the absorption trough 
detected at ~ 9.2 keV is not consistent with the Fexxv K- 
shell edge and the two features are unlikely to be directly 
associated. 

In addition, there are a few other complications which 
could have an effect on the identification of blue-shifted ab- 
sorption profiles at E > 7keV. In particular, because differ- 
ent combinations of xstar parameters (i.e., JVh, log£) can 
yield equivalent solutions to the \ 2 distribution at differing 
red-shifts there can be a level of degeneracy when identi- 
fying the ion responsible for a discrete absorption trough. 
This effect is not particularly significant in sources where 
there are two absorption lines at a common velocity shift 
because the absorber is determined by the joint constraint 
of fitting both profiles, but it can become important in 
the instances where the absorption is manifested through 
a single trough. In these cases it can be difficult to de- 
termine whether an absorption line is due to, for example, 
Fexxv Hea or Fexxvi Lya, which can therefore influence 
the inferred velocity of the outflow and hence any inferred 
absorber energetics. 

There are 10 sources (11 observations) in the sample in 
which a single absorption trough is detected at E > 7 keV at 
the Pmc ^ 95% significance (see Table 2). In each of these 
observations we conducted a search for alternative xstar so- 
lutions by stepping the red-shift of the highly-ionised xstar 
table through the FeK band of each spectrum. This en- 
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abled the x minima for different xstar solutions at differ- 
ent velocity-shifts, for various combinations of column den- 
sity and ionisation parameter, to be mapped. This process is 
analogous to that used by Tf 1 and is useful when it come to 
ascertaining whether there were any alternative fits to the 
absorption lines, and to check for degeneracies between line 
identifications. 

For each source where only a single line is detected we 
took the best-fitting continuum model (including any neces- 
sary soft-band absorbers) and froze all model parameters bar 
those of the highly-ionised absorption table and the normal- 
isation of the primary power-law. To map both blue and red- 
shifted x 2 minima the highly-ionised xstar table was then 
stepped in Az — 10~ 3 increments between —0.5 ^ z < 0.5, 
and contour plots produced after each run. Example con- 
tour plots for Mrk279 and ESO 103-G035, which illustrate 
the cases of non-degenerate and degenerate xstar solutions, 
respectively, are shown in Figure 6. In Mrk279 (Figure 6, 
top panel) the stepping process yields a single valid xstar 
solution at a velocity shift corresponding to Fexxvi Lya. 
However, in ESO 103-G035 (Figure 6, middle panel) there 
are two degenerate x 2 minima found which are statistically 
equivalent at the 90% level (i.e., Ax 2 < 2.71). The mea- 
sured parameters imply that the lowest velocity solution is 
associated with Fexxvi Ka, while the higher velocity solu- 
tion is likely due Fe xxv Hea with some contamination from 
lower ionised species of iron. In this scenario it is difficult to 
unambiguously identify the responsible Fe ion and therefore 
gauge the appropriate outflow velocity. 

Of the 10 sources in which the FeK absorber is man- 
ifested by a single trough only 4 show evidence for degen- 
erate xstar solutions (3C 390.3, 4C +74.26, ESO 103-G035 
and NGC4151). For these cases we report the mean xstar 
parameters as measured from all valid solutions in Table 3 
to account for the uncertainty in line identifications, with 
the associated errors being taken as half of the range be- 
tween the absolute and minimum values. The remaining 7 
sources have only a single x 2 minimum which corresponds 
to the line identification reported in Table 2. In sources with 
degenerate xstar solutions we report the most conservative 
line identification (i.e., associated with Fexxvi Lya) in Ta- 
ble 2. 




-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 
Observed Redshift (zj 

Figure 6. Example plots showing the x 2 statistic versus observed 
red-shift of the xstar absorber (z ). Plots are shown for Mrk279 
(top panel), ESO 103-G035 (middle panel) and NGC5506 (bot- 
tom panel). In all panels the black solid line show the \ 2 con- 
fidence contour and underlying \ 2 distribution. In the bottom 
panel the black and red solid lines show the confidence contours 
for xstar grids both with (dot-dashed red) and without (solid 
black) Nickel included at Solar abundances. In all panels the hor- 
izontal dashed lines show the 90% confidence level for one inter- 
esting parameter, while the vertical dash-dotted line shows the 
best-fitting absorber red-shift. Colour version available online. 



4-6.2 The influence of Nickel 

In principle, the presence of Nickel could also complicate our 
line identifications above 7.1 keV. In particular, the Ka lines 
of Ni xxvil and Nixxvill, which are expected at rest-frame 
energies of E ~ 7.78 keV and E ~ 8.09 keV, respectively, 
could offer a more energetically conservative identification 
for the highest energy absorption lines by virtue of requir- 
ing a lower blue-shifted velocity. A total of six sources have 
absorption lines detected at energies in the regime where 
Nickel could complicate line identifications (Table 2). Be- 
cause of uncertainties in the atomic rates of Nickel 13 its 
abundance is set to zero in xstar by default which means 
that any additional solutions to the x 2 distribution due to 
it would not be picked up during the red-shift stepping. In 
order to robustly search for plausible xstar solutions which 

13 http://heasarc.gsfc.nasa.gov/xstar/docs/html/node74.html 



may be associated with Nickel we therefore re-generated the 
xstar tables using the same assumed parameters as before 
(see Section 4.6), but this time with Nickel included at solar 
abundances, and again searched for alternative xstar solu- 
tions using the method outlined previously. 

The bottom panel of Figure 6 shows a comparison of 
the contour plots obtained for NGC 5506 when Nickel is set 
to zero or Solar in the xstar table. As expected given the 
Z^h/Zfc ~ 0.05 abundance ratio at solar abundances (i.e., 
Grevesse & Sauval 1998), in terms of valid xstar solutions 
for NGC 5506 there are no tangible differences between the 
two grids with Fexxvi Lya being the only valid identifi- 
cation in both cases. Similar is also true for the other 5 
sources with absorption lines detected at E > 7.78 keV, with 
no valid solutions corresponding to Nickel being found in 
any of the cases. Moreover, it should be noted that in order 
to achieve the measured equivalent width of the absorption 
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Table 4. Fraction of detected outflows 



Line(s) 


Sourccs(spectra) 


Fexxv Hea 


2(3) 


Fexxvi Lya 


9(10) 


Fe xxv-xxvi 


7(12) 


Multi u out 


2(3) 


Total 


20 (28) 



lines the required column density would be unreasonably 
large, i.e., log(iVH/cm~ 2 ) > 24, which would absorb the ob- 
served continuum far beyond that which is observed in the 
spectra. 

Therefore, bar the extraordinary case where it is > 20 
times over-abundant relative to solar values, Ni has a negligi- 
ble influence on the Fe K band and identifying the observed 
absorption lines with K-shell resonance lines of Fe is secure. 
We report the best-fitting xstar parameters, including the 
measured outflow velocities relative to the host galaxy, in 
Table 3. The statistical significances as per the A^ 2 and the 
F-test are also reported; in all cases the addition of an xstar 
grid to fit the observed lines improves the fit by at least the 
99% level by the F-test. Note that since some sources have 
either more than one absorption complex detected in a single 
spectrum (i.e., APM 08279+5255 and PDS456), or have ab- 
sorption detected in more than one epoch (i.e., Mrk766), we 
also report the mean absorption parameters on a per source 
basis. The mean values are used in all subsequent analysis to 
prevent sources with multiple FeK absorber detection from 
over- weighting the resultant parameter distributions which 
are discussed in Section 5.2. 



5 RESULTS 

5.1 Line detection rate and phenomenology 

A total of 20/51 sources (in 28/73 fitted spectra; both cor- 
responding to ~ 40% of the total sample) show evidence 
for highly-ionised absorption lines in their Suzaku spectra 
at a Monte Carlo significance of Pmc ^ 95%, with 18/20 of 
these outflows also robustly detected at Pmc ^ 99% signifi- 
cance (see Table 2) . Of the 28 observations with individually 
detected FeK absorbers there are 10 which are consistent 
with having Fexxvi Lya as the dominant Fe ion, 3 where 
Fexxv Hea (and/or lower ionisation species of iron) is the 
main contributor, 12 absorbers with both Fexxv Hea and 
Fe xxvi Lya lines with a common outflow velocity, and a fur- 
ther 3 that have two absorption components with different 
outflow velocities. Taking into account that some sources 
have absorption lines detected in more than one observa- 
tion, on a per-source basis this corresponds to: 2/20 having 
Fexxv Hea absorption, 9 with Fexxvi Lya absorption 14 , 7 
with both Fexxv and Fexxvi Lya absorption, and 2 hav- 
ing multi-Wout systems (out of 20). Interestingly, 8/9 sources 
with just a single Fexxvi Lya absorption line have outflow 

14 In reality the number of Fe xxvi Lya systems could be some- 
what lower than this value due to the 4 AGN where the absorption 
is equally well fitted by Fexxv Hea at a slightly lower velocity 
(i.e., see Section 4.6.1). 



velocities which exceed the v out ~ 10000 km s -1 threshold 
employed by T10A when identifying 'UFO' systems, while 
only 1 of the 6 sources with Fe xxv-xxvi absorption exhibit 
a mean velocity which exceed this threshold value. This is 
consistent with the view that higher ionised outflows origi- 
nate closer to the central AGN and therefore have a higher 
outflow velocity, and is a point which will be further dis- 
cussed in a companion paper (Gofford et al. 2012; in prepa- 
ration). 

Histograms of EW for the absorption lines are shown 
in the top and middle panels of Figure 7, respectively. Note 
that in these plots the individual line profiles in a Fe xxv- 
xxvi pair are considered separately for the sake of clarity, 
and only those lines which have been individually detected 
at a Monte Carlo significance of Pmc ^ 95% are consid- 
ered. The histogram in the bottom panel shows the EW 
distribution for all statistically significant absorption lines, 
including those which are consistent with being a blend of 
the Fe xxv-xxvi Ka transitions. Measured EWs span from 
a few tens of eV up to ~ 130 — 140 eV (Figure 7, bottom 
panel), in a distribution which is consistent with the curve 
of growth analysis conducted by Til. The mean EW of the 
Fexxv Hea and Fexxvi Lya absorption lines are ~ 36 eV 
and ~ 38 eV, respectively, with the total mean EW for all 
detected profiles (i.e., included blended ones) is ~ 40 eV. Im- 
portantly, the observed distribution of absorption line EWs 
is broadly consistent with that found by T10A using XMM- 
Newton (marked with a dashed black line in Figure 7). 

As noted by T10A a useful quantity to calculate is the 
global probability that all of the observed lines are due to 
purely statistical shot noise, which can be done using the 
binomial distribution. For an event with null-probability p 
the chance probability of the event n happening in N trials 
is given by: 

nn-,N, P ) = v ^^ P ^l- P f- (1) 

For n = 20 FeK-band absorption line systems detected in 
N — 51 sources at a Monte Carlo significance of Pmc ^ 95% 
the probability of one of these absorption systems being 
due to shot noise can be taken as p < 0.05. On a per- 
source basis the probability of all of the observed absorption 
systems being associated with noise is then very low, with 
P < 2 x 10~ 13 , which further reduces to P < 5 x 10~ 18 when 
considering the fact that some sources have lines detected 
in more than one epoch. This suggests that the observed 
lines are very unlikely to be associated with simple statis- 
tical fluctuations in the spectra. Moreover, it is important 
to remember that because the vast majority of absorption 
lines are detected at Pmc > 95% significance these probabil- 
ities only represent conservative lower limits on the global 
probability of all lines being false detections. 

5.2 Absorber properties 

The xstar parameter distribution as measured with Suzaku, 
again plotted using the mean absorber parameters averaged 
over all observations, are shown in Figure 8. Overall, the 
general distributions and mean parameter values are broadly 
consistent with those found by T10A using XMM-Newton. 
The Fe K absorbers detected with Suzaku cover a wide range 
of column densities, ranging from 21.5 < log(iY H /cm- 2 ) < 
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Table 3. XSTAR parameters for FeK absorbers. Notes: (1) Source name; (2) Suzaku observation ID; 
(3) Logarithm of absorber column density, in units of cm -2 ; (4) Logarithm of the ionisation parameter, 
in units of erg cm s -1 ; (5) Measured absorber outflow velocity, in units of v/c. Negative outflow veloc- 
ities indicate a net red-shift; (6) Change in x 2 / u when absorber removed from the best-fit model; (7) 
Corresponding absorber significance according to the F-Test. 



Source 


OBSID 


log7V H 


log? 


Vout 


A x /Au 


Pf 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




3C111 


703034010 


23.181°;" 


4 RS+ - 22 


0.039 ±0.005 


18.7/3 


99. 


97 




705040010 


22 -O+0.20 


4 13+ 027 


0.105 ±0.008 


11.6/3 


99. 


08 




(mean) 


99 Qc;+0-28 
zz.yo_g.3g 


A 4^+0.35 
4 -+°-0.23 


0.072 ±0.005 


- 




- 


3C390.3 A 


701060010 


> 23.68 


> 5.46 


0.145 ±0.007 


11.7/2 


> 99. 


11 


4C+74.26 A 


702057010 


> 21.79 


4 06+ 45 


0.185 ±0.026 


14.6/3 


99. 


76 


APM 08279+5255 


stacked [all] 


23.00_°;°£ 


3 4l+°- 08 


0.139 ±0.012 


27.3/3 


> 99. 


99 






23.00* 


3 61+ ' 23 
°- D1 -0.14 


0.431 ± 0.019 


14.3/2 


99. 


59 




(mean) 


23.00_g;il 


4 1-3+023 
4 - ic> -0.17 


0.285 ±0.011 


- 




- 


CBS 126 


705042010 


> 23.73 


4 77 +0.26 
' -0.17 


0.012 ±0.006 


28.9/3 


> 99. 


99 


ESO 103-G035 A 


703031010 


> 21.90 


4 36+ 1 19 
4 - JD -1.19 


0.056 ±0.025 


33.2/3 


> 99. 


99 


MCG -6-30-15 


stacked [all] 


22 16+ 08 
zz. io_ 008 


3 64+ 05 
°-° -0.06 


0.007 ± 0.002 


103.0/3 


> 99. 


99 


MR 2251-178 


705041010 


91 r.A+0-22 


3 26+ 12 


0.137 ± 0.008 


42.4/3 


> 99. 


99 


Mrk 279 


704031010 


23.38_°;« 


4 42+ - 15 
^■^ -0.27 


0.220 ± 0.006 


14.5/3 


99. 


.77 


Mrk 766 


701035010 


22.64l0; 2 2 ? 


4 02+ ' 21 
*- uz -0.13 


0.061 ± 0.008 


14.5/3 


99. 


71 




701035020 


22 76+ ' 10 
zz - ,D -0.13 


q fi7 +0.06 
•3-O<_0.06 


0.017 ± 0.004 


63.2/3 


> 99. 


99 




(mean) 


99 7(1+0-15 
zz - ,u -0.19 


o fi +0.13 


0.039 ±0.006 


- 




- 


NGC 1365 


702047010 


23.92_°;°| 


q aa+0-06 
3-»»_o.o7 


0.014 ± 0.001 


4125.1/3 


> 99. 


99 




705031010 


23.401°;^ 


3.88* 


0.002 ±0.002 


191.8/3 


> 99. 


99 




(mean) 


no 70+0. 06 
zo. (o_ . 06 


q aa+0-06 


0.008 ±0.001 


- 




- 


NGC3227 


703022010 


99 74+ - 11 
zz - '*-0.13 


3 go+0.08 
o -° M -0.11 


< 0.002 


63.1/3 


> 99. 


99 




703022030 


99 c;q+ - 16 
zz. oy_ 20 


q QC+0.13 
•5-80-0.14 


0.007 ± 0.004 


21.0/3 


> 99. 


99 




703022050 


22 62 +0 ' 12 
zz.DZ_ 14 


3 89+° 11 
°-° a -0.12 


0.011 ± 0.004 


37.9/3 


> 99. 


99 




(mean) 


22.66+0 15 


q OO+0.19 
d -°°-0.21 


0.005 ±0.004 


- 




- 


NGC 3516 


100031010 


22.56.0^ 


q O4+O.II 

3.84_ 1Q 


0.004 ± 0.002 


141.9/3 


> 99. 


99 


NGC 3783 


701033010 




3 50+ 013 
• 3 - ou -0.08 


< 0.005 


16.2/3 


99. 


90 




704063010 


91 oq+0.09 
21.S3_o.10 


3 45+ - 06 
o -^ J -0.04 


< 0.008 


48.2/3 


> 99. 


99 




(mean) 


91 7Q+0. 16 
Zl. <y_ 23 




< 0.007 


- 




- 


NGC 4051 


700004010 


99 7Q+0-08 
zz. <o_o.09 


4-05t°° 5 5 


0.020 ± 0.002 


80.3/3 


> 99. 


99 




703023010 


22.80 + o^3 


4 94+0 09 
-0.09 


0.015 ±0.002 


20.8/3 


99. 


99 




(mean) 


22.79l°;i° 


4 00+ - 07 
4 - uu -0.07 


0.018 ±0.001 








NGC4151 A 


701034010 


> 21.74 


3 69+ ' 64 
°- uy -0.64 


0.055 ±0.023 


54.8/3 


> 99. 


99 


NGC 4395 


702001010 


22.841°;^ 


q 2 + - 16 


< 0.001 


19.1/3 


99. 


97 


NGC 5506 


stacked [all] 


9 o 99 +0.19 
ZO.ZZ_ 2g 


5 04+ 29 
°- u -0.17 


0.246 ± 0.006 


16.8/3 


99. 


94 


PDS456 


701056010 


23.04l° ;°l 


4 iq+015 
4 - ly -0.14 


0.253 ±0.008 


20.4/3 


> 99. 


99 






23.04* 


4.19* 


0.292 ± 0.009 


12.8/1 


99. 


68 




705041010 


23.04* 


q Qq + 0.02 

°- ao -0.02 


0.253* 


20.4/3 


> 99. 


99 






23.04* 


3.93* 


0.292* 


11.6/1 


99. 


85 




(mean) 


23 04+ ' 08 
zo.u<±_ 08 


4 06+° 15 
^• UD -0.14 


0.273 ±0.006 








SWJ2127.4+5654 


702122010 


22.78lo:3t 


4 1fi +0 ' 29 


0.231 ± 0.006 


11.5/3 


99. 


24 



* indicates a parameter was tied during spectral fitting; 

A denotes sources with xstar solutions which are degenerate at the 90% level (i.e., the \ 2 statistic for 
the Fexxv Hea and Fexxvi Lya solutions differ by A\ ^ 2.71). In these cases the reported values are 
averaged over the solutions and the errors are inferred as half the range between the absolute maximum 
and minimum values. Significances reported in columns (6) and (7) refer to the least significant of the 
valid solutions. See text for further details. 
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Figure 7. Histogram showing the mean equivalent widths of de- 
tected highly-ionised absorption lines. The top and middle panel 
show the distributions for Fexxv Heo and Fexxvi Lya transi- 
tions, respectively, while the bottom panel shows the total dis- 
tribution of all detected Is— >2p absorption lines including those 
which are detected as a blend of the two Ka transitions. For 
sources with more than one line detected in multiple observa- 
tions the mean equivalent widths have been used. The grey di- 
agonally shaded area and the dashed black line show the distri- 
butions found in this work and with XMM-Newton, respectively. 
The XMM-Newton distributions independently created using the 
values listed in table A. 3 in T10A. 



24.0, with a peak in the distribution at log(7YH/cm -2 ) ~ 
22 — 23. As shown by the dot-dashed (blue) and dotted 
(black) lines, which show the mean log(A^H/cm -2 ) value 
as found with Suzaku and XMM-Newton, respectively, the 
mean is log(A?H/cm _2 ) ~ 23 for both samples. From the 
middle panel the ionisation parameters are in the interval 
2.5 < log(£/erg cm s" 1 ) ^ 6 with the significant fraction of 
Fexxv-xxvi pair systems, which persist over only a rela- 
tively narrow range in ionisation parameter (see Til curve 
of growth analysis), leading to a peak in the distribution 
at log(£/ergcms _1 ) ~ 4. Again, as shown by the verti- 
cal lines the mean ionisation parameter in both samples 
is almost identical at log^/ergcms" 1 ) ~ 4.5. The detec- 
tion of relatively low-ionisation material in the FeK band, 
i.e., log(£/ergcms _1 ) ~ 2.5 — 3.0, is particularly interesting 
and suggests that high velocity absorption could feasibly be 
detected at softer X-ray energies through weak, moderately 
ionised, iron lines. Moreover, the detection of a small fraction 
of absorbers with log(^/ergcms" 1 ) ^ 5 in both this work 
and in the T10A sample raises the possibility that material 




3.5 4 
log(v [1 , II /km s" 1 ) 

Figure 8. Histogram showing the overall distributions of mean 
absorber parameters for each source: (a) logarithm of the mean 
column density; (b) logarithm of the mean ionisation parameter; 
(c) logarithm of the mean outflow velocity. The dot-dashed (blue) 
and dotted (black) vertical lines indicate the mean value of the 
Suzaku and XMM-Newton analyses, respectively. 



may be present in some sources which is so highly ionised 
that even iron is not detectable through spectroscopy. If this 
is the case then the fraction of sources with FeK absorp- 
tion (~ 40%) may represent a lower limit on the number 
of sources with intrinsic nuclear outflows along the line-of- 
sight. 

The log(w out / kms -1 ) distribution (bottom panel) ap- 
pears to be relatively continuous over a broad range of 
velocities; ranging from as low as « out < 1,500 km s _1 
up to twt ~ 100, OOOkms -1 . The median outflow veloc- 
ity is w out ~ 17, OOOkms -1 (~ 0.056c) Ninety per cent 
of the detected outflows have v on t > 1,500 km s' 1 which 
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Table 5. Outflow velocity comparison 



Velocity (kms 1 ) 


Suzaku 


XMM-Newton 


No outflow 


2/20 


2/19 


< Dout C 10, 000 


6/20 


2/19 


i>out > 10, 000 


12/20 


15/19 


v out > 30, 000 


8/20 


9/19 



makes the absorption detected at FeK almost systemati- 
cally faster than the traditional soft X-ray warm absorber. 
Only NGC3227, NGC4395 and NGC3783 have FeK out- 
flows which are consistent with having no outflow velocity. 
From Table 5 the distributions of outflow velocity appear to 
be very similar between both the Suzaku and XMM-Newton 
samples, with the bulk of outflows in both samples having 
v ut > 10,000kms _1 . The Suzaku sample does appear to 
have slightly more low-intermediate velocity systems but, 
given the low number statistics involved, the differences are 
probably not significant. A interesting possibility is that 
both Suzaku and XMM-Newton may be subject to an in- 
strumental bias against the detection of low- velocity absorp- 
tion systems. At the energy resolution of the EPIC and XIS 
CCDs the presence of Fexxv and/or Fexxvi emission lines 
could mask the signature of low-velocity absorption systems 
and thereby introduce a selection effect against their detec- 
tion. The only way to reveal the presence of such absorption 
systems would be through future observations with, for ex- 
ample, the calorimeter aboard Astro-H, which would have 
sufficient resolution to distinguish between individual emis- 
sion and absorption components at low velocity. 

Even so, the similarity between all of the measured 
Suzaku and XMM-Newton parameter distributions can be 
quantitatively assessed using the Kolmogorov-Smirnov two- 
sample test (K-S Test) which uses the maximum differences 
between the cumulative fractional distribution of two sets of 
data to determine the probability that they are both drawn 
from the same parent sample. In this case such a test can 
be used to quantify the level at which the column density, 
ionisation parameter and outflow velocity distributions as 
measured with both Suzaku and XMM-Newton are in agree- 
ment. For a null hypothesis that each of the Suzaku and 
XMM-Newton distributions are drawn from the same par- 
ent sample we are unable to conclusively rule out the null 
hypothesis in any of the three cases at greater than the 90% 
confidence level. 



6 DISCUSSION 

6.1 Detailed sample comparison 

6.1.1 Radio-quiet sources 

In this work we have detected FeK outflows in 17 radio- 
quiet AGN and 3 in radio-loud. Of the detections in radio- 
quiet sources, 9/17 also have had observations included 
in the XMM-Newton outflow sample; with 6/9 of these 
cases having FeK outflows confirmed by the T10A analysis 
(i.e., Mrk766, NGC4051, NGC3516, NGC3783, NGC4151, 
Mrk279). In the remaining 3 sources, namely NGC5506, 
MCG -6-30-15 and NGC3227, the Suzaku outflow detec- 
tions are not confirmed by T10A. Even so, given that in 



the case of NGC 3227 the Fe K absorption lines are detected 
at Pmc > 99.9% significance in 3 separate Suzaku observa- 
tions, and in MCG -6-30-15 the Fe xxv Hea and Fe xxvi Lya 
absorption lines present are both statistically resolved at 
Pmc > 99.9% confidence, the Suzaku line detections alone 
still imply very robust absorber detections in both of these 
sources. A further 3 sources with absorption detected in their 
Suzaku spectra have had statistically significant absorption 
detected in XMM-Newton observations which were not in- 
cluded in the T10A outflow sample (NGC 1365, Risaliti et al. 
2005; PDS456, Reeves et al. 2003, 2009; APM 08279+5255, 
Chartas et al. 2009). Furthermore, weak hard-band absorp- 
tion has been noted in the XMM-Newton spectra of MCG - 
6-30-15 by several authors (e.g., Fabian et al. 2002; Vaughan 
& Fabian 2004; Nandra et al. 2007; Miller et al. 2008). 
In particular, Nandra et al. (2007) reported the presence 
of moderately significant (i.e., just below 99% confidence) 
Fexxv Hea absorption in XMM.OBSID: 00297401010, but 
this detection was not significantly replicated by T10A and 
highlights the need for detailed broad-band spectral models 
when assessing for the presence of highly-ionised absorption 
lines. 

Therefore a total of 9/17 (~ 53%) of the sources with 
outflows detected by Suzaku have also previously been de- 
tected in XMM-Newton data. In addition to these, the ab- 
sorbers in MR 2251-178 and MCG -6-30-15 are also corrobo- 
rated on the basis of their Chandra /HETG data. Gibson 
ct al. (2005) detected a resolved Fexxv Hea absorption 
line at E ~ 7.25 keV in the HETG spectrum MR 2251- 
178 (v out ~ ^OOOkms^ 1 ) while Young et al. (2005) found 
variable absorption lines due to Fexxv and Fexxvi Lycv 
in MCG 6-30-15 with tw ~ 2, 000 km s -1 . In both cases 
the outflow velocities are consistent with those found in 
our Suzaku analysis, and the detection in the HETG spec- 
trum of MCG -6-30-15 overcomes the ambiguities which re- 
mained on the basis of the XMM-Newton data and inde- 
pendently corroborates the Suzaku line detections. FeK ab- 
sorption has also been reported in the Chandra spectrum of 
APM 08279+5255 (Chartas et al. 2002) making it the only 
source in the sample to have its Fe K outflow independently 
detected in all 3 observatories at a high significance level. 

Thus, a total of 11/17 (~ 65%) of the sources 
with FeK outflows detected with Suzaku have also had 
outflows reported elsewhere in the literature with either 
XMM-Newton or Chandra. Of the remaining 6 sources, 
SWIFT J2127.4+5654 has been observed by XMM-Newton 
twice: once in 2009 (~ 24ks) and once in 2010 (~ 131 ks), 
although neither observation has been published at the 
time of writing. CBS 126 have not been observed by either 
XMM-Newton or Chandra, while ESO 103-G035 only has a 
13 ks XMM-Newton observation available which is insuffi- 
cient to test for FeK absorption. This leaves only 3 AGN, 
NGC 5506, NGC 3227 and NGC 4395, with sufficient obser- 
vations available which do not have independent confirma- 
tions with other observatories; although, as discussed previ- 
ously, NGC 3227 does have multiple detections with Suzaku 
which suggests a robust detection. 

6.1.2 Radio-loud sources 

There are 6 Broad-line radio galaxies (BLRGs) in our Suzaku 
sample. This includes all of the sources which were part of 
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the T10B radio-loud UFO case study (i.e., 3C 111, 3C 120, 
3C 445, 3C 390.3 and 3C 382), as well as 4C+74.26 which 
was not part of their analysis. T10B state that 3/5 of 
the BLRGs in their sample, namely 3C 111, 3C 120 and 
3C 390.3, show evidence for highly-ionised absorption lines 
in their Suzaku spectra; with that found in 3C 111 also being 
confirmed in more recent subsequent observations (Tombesi 
et al. 2011) 

We are able to confirm the absorption lines detected 
in both 3C 111 and 3C 390.3 at a similar confidence level. 
However, while we are unable to directly confirm the line 
detection reported in 3C 120 we note that by adding two 
narrow (a = 10 eV) Gaussians at the rest-frame energies for 
the Fexxv Hea and Fexxvi Lya lines reported by T10B 
(7.25 keV and 7.54 keV, respectively) we can place lower 
limits of EWuc a > -8cV and EW Lya > -9eV for the 
Fexxv Hea and Fexxvi Lya absorption lines, respectively. 
These limits are consistent with those reported by T10B at 
the 90% level. Furthermore, we have also detected a high- 
velocity outflow (i>out ~ 0.22 c) in the BLRG 4C+72.26 and 
note that an FeK outflow is also detected in the Chan- 
dra/LETG data for 3C445 (Reeves et al. 2010; see also 
Braito et al. 2011). Thus the number of BLRGs with high- 
velocity outflows could tentatively rise to 4/6 if observations 
from multiple observatories are included, suggesting that 
such outflows could be an important component in radio- 
loud sources. It is important to note, however, that the cur- 
rent sample size of 7 Suzaku observed radio-loud AGN is 
insufficient for a statistical study into the prevalence of Fe K 
absorption in radio-loud sources. That being said, Tombesi 
ct al. (2012; in prep) are currently performing a detailed 
analysis on a large sample of ~ 30 radio-loud AGN observed 
using XMM-Newton and Suzaku which will robustly assess 
for the prevalence of Fe K absorption in these objects. 

6.1.3 Recent Suzaku samples 

Patrick et al. (2012) have recently published the analysis of 
another large sample of Suzaku selected AGN, and is another 
sample with which we can compare our results. While the 
primary aims of the Patrick et al. (2012) study was to assess 
the properties of Fe K emission lines and black hole spin, the 
authors also report the detection of highly- ionised Fe K ab- 
sorption in 14/46 (~ 30%) of their objects. The two samples 
have 43 AGN in common. This includes 15 AGN in which 
we have detected Fe K absorption, of which 11 are confirmed 
by the Patrick et al. (2012) analysis. The four objects where 
our outflow detections are not corroborated are: 3C 390.3, 
4C+74.26, NGC5506 and SWIFT J2127.4+5654. The out- 
flows that we detect in three of these objects (3C 390.3, 
4C+74.26 and SWIFT J2127.4+5654) are only marginal de- 
tections from Monte Carlo simulations, i.e., 95% < Pmc < 
99%, which suggests that the lines may not be visually ap- 
parent in the raw data without a statistically driven analysis. 
The outflow in NGC 5506 is detected with P M c = 99.8% 
which suggests that the discrepancy is likely a result of 
the method in which the datasets were analysed. Indeed, 
our detection of Fexxvi Lya in this source is taken from 
the time-averaged (stacked) Suzaku spectrum while Patrick 
et al. (2012) analysed the observations separately. The out- 
flow in NGC 5506 may be intrinsically weak in the individ- 



ual Suzaku epochs, and only become apparent when they 
are time-averaged due to higher S/N. 

An interesting outcome of the Patrick et al. (2012) anal- 
ysis is that they report outflow detections in three sources 
(3C 445, NGC 5548 and PC 1211+143) which we do not find 
in this work. Outflows in 3C445 and PG 1211 + 143 been de- 
tected with other observatories (3C445: Reeves et al. 2010, 
PG 1211+143: Pounds et al. 2003) and they are known to 
be intrinsically variable, being weakest in their respective 
Suzaku observations (Braito et al. 2011, Reeves et al. 2008b). 
In NGC 5548 we note that while there is a weak absorption 
trough detected at the energy expected for Fe xxv Hea (see 
appropriate entry in Appendix B) it falls below our detec- 
tion criteria. Similar is also true for the Suzaku line detec- 
tion in PG 1211+143, which is only detected at 90% signif- 
icance. However, when compared to the sample of Patrick 
et al. (2012) the broader outcomes are entirely consistent. 
Both analyses find that there is a strong peak in the ionisa- 
tion parameter distribution at log(^/ergcms _1 ) ~ 4 — 4.5, 
with the associated velocities spanning a continuous range 
from ~ < t>out < 100, 000 km s -1 . Therefore, the two sam- 
ples are in complete agreement when it comes to the overall 
properties of the FeK absorbers. 

6.2 Evidence for complex variability and absorber 
structure? 

In addition to broad-band spectral variability there is also 
compelling evidence for intricate variability in the FeK 
absorber itself. In a number of cases this variability can 
be simultaneously and self-consistently fitted using xstar 
absorption tables which strongly implies an atomic origin 
rather than an association with spectral fluctuations and 
photon statistics. In NGC 3227 transient FeK absorption 
lines are detected in 3 of the 6 Suzaku epochs, with the line 
parameters varying on a At ~ 2 week time-scale. Similar 
is also true in Mrk766, PDS456 and NGC 1365, which all 
exhibit absorption line variability albeit over a longer (i.e., 
At ~year) baseline. In both Mrk766 and PDS456 the ab- 
sorber appears to change in ionisation states between the 
Suzaku epochs, with the outflow in Mrk766 also showing a 
decrease in outflow velocity. Here, we discuss both sources 
in greater detail. 

Mrk766 and PDS456 have both long been known to 
harbour Fe K absorbers (Mrk 766: Pounds et al. 2003; Turner 
et al. 2007; Miller et al. 2008; Risaliti et al. 2011. PDS456: 
Reeves et al. 2000, 2003, 2009). The range of rest-frame en- 
ergies for the absorption profiles measured with Suzaku in 
each source are consistent with those previously published in 
the literature. In Mrk 766 we find a common velocity shift 
of v out ~ 6,000 - 18,000kms _1 for the Fexxv Hea and 
Fexxvi Lya lines, which is consistent with the iw.xmm ~ 
3, 000-16, 000 km s" 1 reported previously in XMM-Newton 
data (Miller et al. 2007; Turner et al. 2007; Risaliti et al. 
2011). Similarly, from our simultaneous analysis of the 2007 
and unpublished 2011 Suzaku observations of PDS456 we 
find two absorption troughs which persist across both se- 
quences (see Figure 9). These troughs - which are likely due 
to blends of the Fe xxv Hea and Fe xxvi Lya lines owing to 
their breadth - have a mean inferred {v out } ~ 0.3 c, which 
is consistent with that reported by Reeves et al. (2009) in 
their initial analysis of the 2007 observation. 
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An interesting property of the outflows in both of these 
sources is that the parameters of the hard-band absorber 
appears to be linked to both the incident source flux, and 
to the broad-band absorption characteristics of the X-ray 
spectrum. On the basis of detailed time- and flux-resolved 
spectroscopy of Mrk766 using XMM-Newton, Turner et al. 
(2007) noted that the EW of the Fe xxvi Lya absorption 
line could be correlated with the intrinsic source flux, with 
the correlation not easily accounted for by simply appeal- 
ing to changes in the ionisation state of a constant column 
density absorber in response to changes in the incident con- 
tinuum. This provided evidence for complex absorber struc- 
ture in Mrk766. Further evidence for Fe K absorber com- 
plexity was suggested by Risaliti et al. (2011) who, from a 
re-analysis of the same XMM-Newton data, found that the 
presence of the highly-ionised absorption profiles appeared 
to be contingent on the source being eclipsed by a number 
of high column density, i.e., log(A r H/cm -2 ) > 23, partially- 
covering clouds with a low ionisation state. This led the 
authors to propose a scenario where multiple 'cometary' ab- 
sorption clouds with a large gradient in both column density 
and ionisation state were moving across the line of sight, thus 
giving rise to the varying absorption characteristics present 
in different time-slices (e.g., v out ). 

Our Suzaku analyses yields further evidence for a com- 
plex and dynamic structure in hard-band absorbers. In 
Mrk 766 the EW of the Fe xxv Hea and Fe xxvi Lya lines 
present in the Suzaku spectrum are statistically consistent 
at the 90% level but, as shown in the appropriate entries in 
Appendix B, there is a visual shift in the dominant iron ion 
of the absorber, with Fexxvi Lya (EWvexxvi = — 40±17eV) 
appearing strongest in OBSID 701035010 and Fexxv Hea 
(EW Fcxxv = -60 ± 14 eV) strongest in OBSID 701035020, 
which is possibly associated with the A log(£/ergcm s _1 ) ~ 
0.4 change in ionisation parameter between the two obser- 
vations. Furthermore, and as mentioned previously, the de- 
crease in ionisation is also met with a decrease in measured 
outflow velocity and the onset of strong Fexxv emission 
which could be associated with a cloud/clump of material 
moving outside of the sight-line. A similar effect is also 
apparent in PDS 456 where, despite there being no mea- 
surable change in the relative velocity of the two blended 
Fexxv-xxvi absorption complexes, the absorber ionisation 
state decreases by A log(£/ergcm s _1 ) ~ 0.26 in the 2011 
observation relative to 2007 which leads to a slight visual 
broadening of the absorption profiles due to an increased 
Fe xxv Hea contribution in the blended absorption line pro- 
file. Interestingly, the decreased ionisation state of the ab- 
sorbers in both Mrk 766 and PDS 456 occurs when the cov- 
ering fraction of the soft X-ray absorber is at its highest 
which implies a dynamic clumpy structure to the outflow. 
Moreover, as noted by both Turner et al. (2007) and Behar 
et al. (2010), respectively, the opposite is also true in both 
sources, i.e., when the sources are in their least absorbed 
state, the discrete highly-ionised absorption lines are not 
significantly present in the X-ray spectra of either source. 
This is consistent with the conclusions reached by Risaliti 
et al. (2011) in the context of stratified, partially-covering, 
absorption clumps needing to be present along the line of 
sight, and close to the continuum source, for highly-ionised 
absorption lines to be observed. 

Additional evidence for this supposition is provided 
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Figure 9. Ratio plot comparison of the 2007 (black squares) and 
2011 (red circles) Suzaku observations of PDS 456 in the FeK 
band. The dashed vertical lines indicate the position of the two 
blended Fe xxv-xxvi absorption troughs present in both spectra, 
at rest-frame energies of ~ 9.0 and ~ 9.6 keV, respectively. There 
is a noticeable broadening in the left hand absorption profile in 
the 2011 observation which is likely caused by a drop in absorber 
ionisation state, and an associated increase in the Fexxv Hea 
contribution to the absorption profile. See text and Reeves et al. 
(2012, in prep) for further details. Colour version available online. 



through consideration of our spectral models to other 
sources in the sample. It is intriguing to note that statis- 
tically acceptable (i.e., Xreduccd ~ 1-0) partially covering ab- 
sorption models are obtained for a curiously large fraction 
of the sources with detected outflows (12/20; ~ 60%) which 
further suggests that there may be a link between complex 
absorption geometries and highly-ionised absorption lines. 
One possibility is where the partially covering components 
represent denser clumps of material in an inhomogeneous 
highly-ionised wind from the accretion disc, similar to the 
case reported in Mrk 766 (Miller et al. 2007; Turner et al. 
2007; Risaliti et al. 2011). Indeed, inhomogeneous winds 
with stratified or filamentary ionisation and density struc- 
ture are expected as a natural consequence of the accretion 
process and are ubiquitously seen in both hydrodynamic 
(e.g., Proga & Kallman 2004; Kurosawa & Proga 2009; Sim 
et al. 2010) and magneto-hydrodynamic (e.g., Ohsuga et al. 
2009; Ohsuga & Mineshige 2011) simulations of accretion 
discs. It is therefore not particularly surprising to find a 
possible link between clumpy absorption dominated models 
and sources with high- velocity outflows. 

6.3 On the claimed publication bias 

The transient and variable nature of FeK absorption lines 
has in part lead to their true veracity to be questioned in 
the literature. In particular, Vaughan & Uttley (2008) sug- 
gested that there may be a publication bias at play in the 
reporting of both red- and blue-shifted features in the Fe K 
band, with only those observations with the strongest line 
detections being reported in the literature. Through a plot 
of '£H/-error(EW)' Vaughan & Uttley (2008) showed that 
the £W/error(EW) ratio remained relatively constant (i.e., 
lines with larger EW have correspondingly larger 90% er- 
rors) over a wide range of EWs in their sample of narrow 
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Figure 10. Plot of the measured equivalent width of the blue- 
shifted absorption lines with respect to their negative 90% error. 
The circles (red), diamonds (green), and squares (blue) corre- 
spond to the Fe xxv Hea, Fe xxvi Lya and blended Fe xxv-xxvi 
lines, respectively. Analogous values for the blue-shifted absorp- 
tion lines found in the T10A XMM-Newton sample are shown as 
pentagons (grey). The diagonal dashed line shows the expected 
trend if EW = — error(EW) and the dotted (dark-cyan) area 
shows the area of non-detection (see text for further details). 
Colour version available online. 

lines in 38 sources collected from the literature. The au- 
thors also noted that since the vast majority of reported 
detections (prior to 2008, at least) were of relatively low sig- 
nificance (i.e., typically around 2 — 3a) suggested that the 
lines may be spurious, and more consistent with merely be- 
ing the strongest natural fluctuations in otherwise feature- 
less spectra than real atomic features. Moreover, Vaughan & 
Uttley (2008) suggested that the conspicuous absence of any 
lines in the upper left quadrant of the EW — error (EW) plot 
(which would correspond to stronger lines having smaller un- 
certainties) implied that the detection of the narrow velocity 
shifted lines was in some way inversely correlated to the sta- 
tistical quality of the observation. Observations with longer 
exposures and better photon statistics only ever showed the 
weakest lines, which further enforced the possibility of a bias 
in the published observations. 

However, these points were addressed through the sys- 
tematic search conducted by T10A using XMM-Newton, at 
least in the context of blue-shifted FeK absorption lines. 
By uniformly and comprehensively searching for such fea- 
tures in a complete sample of AGN observations, carefully 
reporting the fraction of detections to null detections, and 
assessing the statistical significance of any detected absorp- 
tion lines through Monte Carlo simulations, T10A were able 
to robustly assess the fraction of AGN in the local uni- 
verse that have Fe K outflows in a way which overcame any 
publication biases and accounted for the possibility of ran- 
dom fluctuations in the source spectra. Our work with the 
Suzaku outflow sample complements the findings of T10A, 



Figure 11. Plot of line significance, <r (= EW/— error(EW)), 
versus total counts between 5 — lOkeV. The blue circles, red 
squares and green crosses show the positions of Fexxv Hea, 
Fexxvi Lya and blended Fexxv-xxvi lines, respectively. Only 
absorption lines detected at Pmc ^ 95% are shown. The dashed 
grey band shows the 2 — 3<r detection regime. Colour version avail- 
able online. 

and lends further weight to the assertion that such out- 
flows are an important intrinsic feature of the AGN X- 
ray spectrum in a large fraction of sources. In Figure 10 
we replicate the Vaughan & Uttley (2008) £W/error(EW) 
plot 15 used by T10A (see their Figure 8) to include data 
points corresponding to the absorption lines found in this 
work at a Pmc ^ 95% significance level. The figure con- 
cisely shows that the distribution of points from both the 
Suzaku and XMM-Newton samples diverges from that of 
the EW=error(_BM / ) 'detection line', with both sets of data 
having a number of points which veer towards the important 
upper left quadrant of the diagram indicating stronger lines 
with smaller uncertainties. Interestingly, the overall distri- 
bution of points for the Suzaku lines follows a similar trend 
to those obtained with XMM-Newton which implies that 
all points are drawn from the same parent population, as 
would be expected should we be studying the same physical 
phenomenon which imprints real spectroscopic lines in the 
X-ray spectrum. 

A further means of visually assessing for a publication 
bias, and for testing whether the absorption lines are con- 
sistent with random fluctuations or real features, is through 
a plot of line significance versus photon counts in the Fe K 
band (i.e., between 5 — lOkeV). Such a plot is shown in Fig- 
ure 11. For simplicity we grade line significance according to 
Na = EW/— error(EW), where it is important to recognise 

15 Note that the axes in this plot use a linear scale rather than the 
logarithmic one as used by Vaughan & Uttley (2008). As noted by 
T10A the use of a logarithmic scale visually compresses the data 
points towards the EW=error(i?lT) threshold of non-detection. 
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that in this case — error(EW) refers to the la (68%) nega- 
tive errors on an absorption lines EW, which were measured 
independently for each absorption line, and not to the 90% 
errors (w 1.6a) reported in Table 2. By using the la negative 
errors the standard Gaussian probabilities (i.e., 2a = 95.5%, 
3a = 99.73%, 4a = 99.994%) then correspond to the statisti- 
cal significance of a line relative to it having an intrinsic EW 
of OeV, and thus not being detected. Figure 11 shows that 
the detected lines typically have significances ranging from 
~ 2a through to ~ 10a over a wide range in count rate; with 
a significance of ~ 30a the extremely strong Fe xxv Hea line 
in NGC 1365 (OBSID 702047010) is by far the most signif- 
icant line in the sample and is omitted from the main plot 
for clarity. 

In such a plot we would expect to see a correlation be- 
tween total Fe K counts and the significance of an absorption 
line, where observations with higher counts in the Fe K band, 
and therefore better photon statistics, have stronger absorp- 
tion line detections. Indeed, while there is noticeable scatter 
in Accounts versus a - as would be expected given that there is 
a wide dispersion in absorption properties (i.e., jV H , £, v out ) 
for a given total counts - a simple Spearmans rank associ- 
ation shows that the data are positively correlated (Spear- 
mans rank coefficient, r s — 0.4406). This is sufficient to rule 
out the possibility of the two parameters being completely 
independent at the > 99% significance for 34 matched pairs. 
Therefore, a values tend to increase in higher counts obser- 
vations which causes the data to gradually diverge from the 
2 — 3a significance region towards the upper right quadrant 
of the plot. Moreover, the vast majority of the line detec- 
tions are located well away from the 2 — 3a 'noise' level, as 
would be expected if the majority of the line detections are 
not purely due to photon statistics. 



7 SUMMARY & CONCLUSIONS 

Making use of data from the Suzaku Data Archives and 
Transmission System (DARTS) we have constructed broad- 
band spectral models for 51 Type 1.0-1.9 AGN to search for 
the presence of Fexxv Hea and/or Fexxvi Lya absorption 
lines at E > 6.7 keV, robustly assessed the statistical sig- 
nificance of any detected absorption systems using detailed 
Monte Carlo simulations, and probed the properties of the 
absorbing material using the xstar (v. 2.2. lbc) photoionisa- 
tion code. The primary findings of this work are as follows: 

1. We find that 20/51 sources (in 28/73 fitted spectra) 
have statistically significant Fexxv Hea or Fexxvi Lya 
absorption lines at E ^ 6.7 keV in the source rest-frame, 
which corresponds to w 40% of the sample. 18 of the 
20 FeK absorption complexes are robustly detected at 
Pmc 99% from Monte Carlo simulations, with those 
remaining only narrowly falling short of the 99% criterion 
required for a robust detection. All absorption lines are 
detected independently and simultaneously in two (or 
more) of the available XIS detectors which further en- 
forces the robustness of the line detections as real features 

2. The detected outflows fall into a range of phenomenolog- 
ical categories. Absorption due to a pair of Fexxv Hea 
and Fexxvi Lya from the same photoionised absorber 
account for 7/20 of the detected outflows (in 12/28 



spectra) making them the most frequently observed form 
of highly-ionised outflow in the sample. Absorption due 
to solely Fexxv Hea or Fexxvi Lya is less frequently 
observed, being detected in 2 sources (3 spectra) and 
9 sources (10 spectra), respectively. The remaining 
two sources have outflows which are best fitted by two 
unresolved and blended Fe xxv-xxvi absorption systems 
with different outflow velocities. 

3. By fitting the absorption with the xstar photoion- 
isation code we find that the absorbers are char- 
acterised by 21.5 < log(A?H/cm~ 2 ) < 24 and 
2.5 < log(^/ergcms _1 ) ^ 6, with mean values of 
log(AfH/cm~ 2 ) « 23 and log(^/ergcms _1 ) « 4.5, re- 
spectively. The distribution of outflow velocities cover 
a wide range, from v ou t < l,500kms _1 up to v out ~ 
100, OOOkms -1 , with 90% of the observed absorbers hav- 
ing v out > 1,500 km s _1 which makes them systemati- 
cally faster than typical absorbers found in the soft X- 
ray band. The median outflow velocity is of the order 
of ~ 0.056 c. Moreover, 60% of the absorbers also have 
Vout > 10, 000 km s -1 and are thus consistent with the 
so called 'Ultra Fast Outflows' (UFOs) defined by T10A. 
A Kolmogorov-Smirnov two-sample test shows that the 
overall distributions for all three parameters are not sta- 
tistically distinguishable from the analogous distributions 
measured by T10A and Til using XMM-Newton. 

Overall, our work with the Suzaku sample is consistent with 
those obtained using XMM-Newton and, combined, the two 
studies provide strong evidence for the presence of very 
highly-ionised, often high-velocity, outflows in the central re- 
gion of a large fraction of not only radio-quiet AGN, but also 
suggests that they may be prevalent in radio-loud sources as 
well. The possible prevalence of very highly-ionised, high- 
velocity winds in AGN is consistent with theoretical mod- 
els which argue that such outflows are an important phe- 
nomenon which may play a role in Galactic-scale feedback 
scenarios. We return to this point in a forthcoming work 
(Paper II; Gofford et al., in prep) where we probe the global 
energetics of the detected absorbers, assess for any correla- 
tions between the outflow parameters and those of the cen- 
tral SMBH, and discuss both the origins of the absorbing 
material and its likely launching mechanism. 
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APPENDIX A: SUMMARY OF OBSERVATIONS 

Table Al. Summary of observation parameters. Notes: (1) Source name; (2) Source type; (3) Cosmological red-shift according 
to the NED; (4) Galactic absorption column density, in units of xl0 21 cm~ 2 ; (5) Suzaku observation ID; (6) Starting date of 
observation; (7) XIS(HXD) net exposures, in kilo-seconds; (8) Total number of (background subtracted) 2 — lOkcV counts for 
the (co-added) XIS-FI(BI) spectra, in units of xlO 3 ; (9) Total background subtracted PIN counts in the 15 — 50keV band, in 
units of XlO 3 ; 



Source 


Type 


red-shift 




OBSID 


Date 


Exposure 


Source net counts 






(z) 




yyyy-mm-dd 


XIS(PIN) 


XIS-FI(BI) 


PIN 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1H 0419-577 


Sy 1.5 


0.10400 


0.126 


702041010 


2007-07-25 


205.9(142.6) 


91 1 4(1 1 1 8 s ! 


49.8 










704064010 


2010-01-16 


246.2(122.8) 


SQ fifAf. x.) 

O i7. \J l rt(J . 'J I 


36.6 










stackcd[all] 






302 9(4 60 41 


79.6 


3C 111 


BLRG 


0.04850 


2.910 


703034010[a] 


2008-08-22 


122.4(101.9) 


106.3(49.8) 


36.2 










705040010[b] 


2010-09-02 


80.7(67.9) 


199 7(95 41 

ICC # 1 1 Jtjt^ 1 


26.6 










705040020 [c] 


2010-09-09 


79.4(66.5) 


258.8(127.1) 


27.4 










705040030 [d] 


2010-09-16 


80.4(65.6) 


9^1 0(1 97 9 s ! 


28.1 


3C 120 


BLRG 


0.03300 


1.060 


700001010[a]t 


2006-02-09 


41.9(31.9) 


1 52 3(50 3") 


13.1 










700001020[b]t 


2006-02-16 


41.6(34.5) 


1 1A ^(44 7~\ 

104.0141.1 J 


15.5 










700001030[c]t 


2006-02-23 


40.9(36.2) 


1 1(4-^ fil 

±- <J kj . 'J y-tij . I 


15.7 










700001040[d]t 


2006-03-02 


40.9(37.9) 


loy.o^io.o j 


16.1 










stacked [bed] 






WC. 11 lOO.O J 


44.6 


3C 382 


BLRG 


0.05890 


0.698 


702125010 


2007-04-27 


130.6(114.3) 


274 5(1 36 1 1 

I I.OI IOU.1 1 


45.5 


3C 390.3 


BLRG 


0.05610 


0.347 


701060010 


2006-12-14 


99.4(92.1) 


151 6(71 81 


38.3 


3C445 


BLRG 


0.05590 


0.449 


702056010 


2007-05-25 


139.8(109.5) 


41 9(20 21 


40.8 


4C +74.26 


BLRG 


0.10400 


1.160 


702057010 


2007-10-28 


91.6(87.3) 


107 0(36 51 


33.8 


APM 08279+5255° 


QSO 


3.91000 


0.411 


701057010[a]* 


2006-10-12 


102. 3(-) 


4.6(1.6)° 












701057020[b]* 


2006-11-01 


102. 3(-) 


4.5(1.6)° 












701057030[c] 


2007-03-24 


117.2R 


3.3(1.8)° 












stacked [all] 










Ark 120 


Sy 1.0 


0.03270 


0.978 


702014010 


2007-04-01 


100.9(89.5) 


1 ^9 1 (7R 7\ 

1 OA .11 1 O. 1 J 


32.3 


Ark 564 


NLSyl 


0.02468 


0.564 


702117010 


2007-06-26 


100.0(81.3) 


1 09 7(*\A A\ 


23.5 


CBS 126 


Sy 1.2 


0.07888 


0.097 


705042010 


2010-10-18 


101.5(84.2) 


91 fiMn 7~\ 

Zr 1 .U ^ 1U . 1 J 


24.8 


ESO 103-G035 


Sy 1.9 


0.01329 


0.764 


703031010 


2008-10-22 


91.5(75.1) 


77 Sl"31 


34.2 


Fairall 9 


Sy 1.2 


0.04700 


0.316 


702043010 


2007-06-07 


167.8(127.3) 


227.3(117.7) 


49.3 










705063010 


2010-05-19 


229.3(162.2) 


226 9(1 03 91 


56.8 


IC4329A 


Sy 1.2 


0.01610 


0.461 


702113010[a] 


2007-08-01 


25.5(20.1) 


131.0(64.8) 


11.4 










702113020[b] 


2007-08-06 


30.6(24.1) 


1 Q7 8(1 00 1 1 
ly i .o^iuu.i j 


15.4 










703113030[c] 


2007-08-11 


24.2(22.1) 




14.5 










703113040[d] 


2007-08-16 


24.2(18.9) 


1 11 1(f\7 1 1 

lOO.O^U '■-'-/ 


11.5 










703113050[e] 


2007-08-20 


24.0(17.5) 


R1 9(40 


9.4 










stackcd[all] 






71 9 7(1^Q (W 


61.7 


IGR J16185-5928 


NLSyl 


0.03500 


2.070 


702123010 


2008-02-09 


76.7(69.6) 


24.1(8.1) 


15.9 


IGR J21247+5058 


Sv 1 


0.0200 


10.000 


702027010 


2007-04-16 


85.0(66.4) 


911 1 (1 Cil 7 s ! 

^OO.l^lUO. 1 J 


38.4 


MCG -02-14-009 


Sy 1.0 


0.02845 


0.948 


703060010 


2008-08-28 


142.2(120.0) 


24.1(10.6) 


20.5 


MCG -2-58-22 


Sy 1.5 


0.04686 


0.291 


704032010 


2009-12-02 


139.0(98.0) 


324.0(129.6) 


38.3 


MCG -5-23-16 


Sy 1.9 


0.00849 


0.800 


700002010 


2005-12-07 


95.7(79.7) 


626.8(198.5) 


54.0 


MCG -6-30-15 


Sy 1.2 


0.00775 


0.392 


700007010[a]* 


2006-01-09 


143.3(118.9) 


504.2(167.0) 


51.3 










700007020[b]t 


2006-01-23 


98.5(76.8) 


304.6(105.0) 


35.1 










700007030[c]* 


2006-01-27 


96.7(89.8) 


327.5(109.1) 


38.7 










stacked [all] 






1136.3(381.1) 


120.3 


MCG +8-11-11 


Sy 1.5 


0.02050 


1.840 


702112010 


2007-09-17 


98.8(82.9) 


295.8(141.2.2) 


42.1 


MR 2251-178 


RQQ 


0.06400 


0.242 


704055010 


2009-05-07 


136.9(103.8) 


330.9(171.6) 


38.7 


Mrk 79 


Sy 1.2 


0.02220 


0.527 


702044010 


2007-04-03 


83.7(76.9) 


70.5(36.6) 


20.1 


MrkllO 


Sy 1.0 


0.03530 


0.130 


702124010 


2007-11-02 


90.9(80.4) 


93.7(46.4) 


29.4 


Mrk 205 


Sy 1.0 


0.07085 


0.240 


705062010 


2010-05-22 


101.5(85.3) 


36.0(15.9) 


29.3 


Mrk 279 


Sy 1.0 


0.03050 


0.152 


704031010 


2009-05-14 


160.4(139.8) 


43.8(23.5) 


40.3 


Mrk 335 


NLSyl 


0.02580 


0.366 


701031010* 


2006-06-21 


151.3(131.7) 


162.6(49.0) 


45.9 


Mrk 359 


Sy 1.5 


0.01739 


0.426 


701082010 


2007-02-06 


107.5(96.1) 


27.2(13.6) 


8.8 



*Only the XIS 3 spectrum was available during the observation of SWIFT J2127. 4+5654; 

' Indicates that the XIS 2 was operational during the observation and is included as part of the co-added XIS-FI spectrum.; 
a The total net source counts listed in column (8) for APM 08279+5255 are taken from 0.6 — 10 kcV in the observer frame. See 
text for further details. 
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Table Al - continued - Summary of observation parameters. 



Source 


Type 


red-shift 


jyGal 


OBSID 


Date 


Exposure 


Total counts 






(*) 






yyyy-mm-dd 


XIS(PIN) 


XIS-FI(BI) 


PIN 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(°) 


(9J 


Mrk 509 


Sy 1.5 


0.03440 


0.425 


701093010[a]t 


2006-04-25 


24.6(14.5) 


93.8(30.8) 


6.6 










701093020 [b]* 


2006-10-14 


25.9(21.2) 


122.9(43.8) 


10.5 










701093030 [c] 


2006-11-15 


24.5(17.3) 


73.9(40.8) 


8.0 










701093040 [d] 


2006-11-27 


33.1(27.6) 


92.4(49.3) 


12.5 










stacked [bed] 


- 


- 


289.3(133.8) 


30.8 










705025010 


2010-11-21 


102.1(85.7) 


325.2(171.4) 


37.8 


Mrk 766 


Sy 1.0 


0.01290 


0.178 


701035010 


2006-11-16 


97.9(90.5) 


175.5(108.5) 


31.6 










701035020 


2007-11-17 


59.4(47.7) 


73.0(40.3) 


17.1 


Mrk 841 


Sy 1.5 


0.03640 


0.222 


701084010[a] 


2007-01-22 


51.8(43.7) 


36.3(18.0) 


14.5 










701084020 [b] 


2007-07-23 


50.9(44.4) 


35.9(18.4) 


16.1 










stacked [all] 


— 


— 


72.2(37.4) 


31.3 


NGC1365 


Sy 1.8 


0.00546 


0.134 


702047010 


2008-01-21 


160.5(136.6) 


96.8(48.2) 


53.9 










705031010 


2010-06-27 


151.6(114.3) 


35.5(15.9) 


42.2 










705031020 


2010-07-15 


302.2(231.5) 


45.8(21.4) 


76.2 


NGC 2992 


Sy 1.9 


0.00771 


0.487 


700005010[a]* 


2005-11-06 


38.8(29.9) 


23.0(6.7) 


9.8 










700005020 [b]* 


2005-11-19 


39.7(31.9) 


39.6(13.1) 


10.4 










700005030 [c]* 


2005-12-13 


46.9(41.5) 


41.9(12.8) 


14.4 










stacked [all] 


- 


- 


84.7(25.0) 


36.1 


NGC3227 


Sy 1.5 


0.00386 


0.199 


703022010 


2008-10-28 


58.9(48.1) 


117.7/60.8 


22.1 










703022020 


2008-11-04 


53.7(46.7) 


38.5(17.3) 


19.6 










703022030 


2008-11-12 


56.6(46.7) 


61.0(28.5) 


20.4 










703022040 


2008-11-20 


64.6(43.4) 


24.1(11.8) 


15.1 










703022050 


2008-11-27 


79.4(37.4) 


69.7(32.9) 


15.4 










703022060 


2008-12-02 


51.4(36.9) 


32.1(14.7) 


15.4 


NGC 3516 


Sy 1.5 


0.00884 


0.345 


100031010* 


2005-10-12 


134.6(115.4) 


233.3(75.4) 


50.7 










704062010 


2009-10-28 


251.4(178.2) 


149.0(73.4) 


55.6 


NGC 3783 


Sy 1.5 


0.00973 


0.991 


701033010* 


2006-06-24 


75.7(68.3) 


302.4(103.0) 


32.2 










704063010 


2009-07-10 


172.8(174.0) 


472.1(234.0) 


84.0 


NGC 4051 


NLSyl 


0.00234 


0.115 


700004010* 


2005-11-10 


119.6(112.6) 


76.4(24.3) 


42.0 










703023010 


2008-11-06 


274.5(204.5) 


349.6(132.5) 


72.0 










703023020 


2008-11-23 


78.4(58.5) 


72.0(35.0) 


21.3 


NGC 4151 


Sy 1.5 


0.03320 


0.230 


701034010 


2006-12-18 


125.0(123.5) 


218.2(97.5) 


75.8 


NGC 4395 


Sy 1.8 


0.00106 


0.135 


702001010 


2007-06-02 


101.5(95.0) 


19.4(9.2) 


30.4 


NGC 4593 


Sy 1.0 


0.00900 


0.189 


702040010 


2007-12-15 


118.8(101.6) 


66.3(25.1) 


34.7 


NGC 5506 


Sy 1.9 


0.00618 


0.408 


701030010[a]* 


2006-08-08 


47.8(38.6) 


372.5(119.7) 


25.4 










701030020 [b]* 


2006-08-11 


53.3(44.8) 


438.9(141.2) 


29.9 










701030030 [c] 


2007-01-31 


57.4(44.7) 


281.1(133.5) 


29.1 










stacked [all] 


— 


— 


1067.5(392.4) 


82.2 


NGC 5548 


Sy 1.5 


0.01718 


0.155 


702042010[a] 


2007-06-18 


31.1(25.6) 


11.0(5.3) 


7.5 










702042020 [b] 


2007-06-24 


35.9(31.2) 


21.7(10.8) 


9.3 










702042040 [c] 


2007-07-08 


30.7(27.0) 


35.9(16.4) 


8.5 










702042050 [d] 


2007-07-15 


30.0(24.5) 


23.3(10.7) 


7.7 










702042060 [e] 


2007-07-22 


28.9(23.1) 


43.3(19.9) 


7.9 










702042070 [f] 


2007-07-29 


31.8(27.6) 


31.6(14.5) 


8.5 










702042080 [g] 


2007-08-05 


38.8(30.4) 


20.5(10.0) 


9.6 










stacked [all] 


— 


— 


187.4(107.8) 


66.3 


NGC 7213 


Sy 1.0 


0.00584 


0.106 


701029010* 


2006-10-22 


90.7(84.3) 


370.2(145.0) 


35.3 


NGC 7469 


Sy 1.0 


0.01630 


0.445 


703028010 


2008-06-24 


112.1(85.3) 


116.2(58.2) 


32.3 


PBCJ0839. 7-1214 


QSO 


0.19787 


0.568 


705007010 


2010-05-08 


80.6(-) 


35.5(17.7) 


_ 


PDS 456 


RQQ 


0.18400 


1.960 


701056010 


2007-02-24 


190.6(-) 


34.0(16.7) 












705041010 


2011-03-16 


125.5(-) 


17.0(7.9) 




PG 1211+143 


RQQ 


0.08090 


0.274 


700009010* 


2005-11-24 


97.6(78.8) 


19.4(5.9) 


- 


PKS 0558-504 


NLSyl 


0.13720 


0.337 


701011010[a] 


2007-01-17 


20.6(17.6) 


11.8(5.9) 


4.8 










701011020 [b] 


2007-01-18 


18.9(16.7) 


16.0(8.1) 


6.5 










701011030[c] 


2007-01-19 


21.3(17.8) 


12.0(6.2) 


6.8 










701011040[d] 


2007-01-20 


19.7(16.0) 


18.3(8.2) 


6.3 










701011050[e] 


2007-01-21 


19.5(15.2) 


17.1(8.2) 


5.4 










stacked [all] 






76.2(44.6) 


21.0 


RBS 1124 


Sy 1.2 


0.20800 


0.152 


702114010 


2007-04-14 


86.2(83.0) 


16.6(7.4) 


21.5 


SWJ2127.4+5654 


Sy 1.0 


0.01440 


7.650 


702122010* 


2007-12-09 


91.7(83.3) 


69.5(81.4) 


31.2 


TONS180 


NLSyl 


0.06198 


0.136 


701021010 


2006-12-09 


120.7(102.4) 


43.5(-) 


27.2 
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APPENDIX B: RATIO AND CONTOUR PLOTS 



H 0419-577 (OBSID:7020410 



3C1 1 1 (OBSID: 703034010) 



3C1 1 1 (OBSID: 7 



Iff 




; i 

III III 1 






ffr 











ime Energy (keV) 
BSID: 705040030) 



Rest-frame Energy (keV) 
3C390.3 (OBSID: 701060010) 



Rest-frame Energy (keV) 
3C445 (OBSID: 70205601 0) 



Rest-frame Energy (keV) 
APM0827S+5255 (OBSID: 7010570[1 0+20+30]) 



II 

Iff ; !';if 





1 










J ^JL 







sl-frame Energy (keV) 
4 (OBSID: 702117010) 



Rest-frame Energy (keV) 
CBS126 (OBSID: 705042010) 



Resl-frame Energy (keV) 
IGRJ21 247+5058 (OBSID: 70202701 0) 





j| i 1 


||tp|j f 




^§4 







■i4 1 






















5 6 


est-trame Energy (keV) 


MCG-2 


-14-9 (OBSID: 703060010) 



hSC 


7 8 9 
Rest-frame Energy (keV) 

103-G035 (OBSID: 703031010) 


1 

illl 1 ! 


| 

Iff 




pi! f 










L_ - . 

















'ML * 




6 


Resl-frame Energy (keV) 


IGRJ16 


185-5928 (OBSID: 702123010) 



Rest-frame Energy (keV) 
MCG-5-23-1 6 (OBSID: 700002010) 



Rest-frame Energy (keV) 
MR 2251-178 (OBSID: 704055010) 





Figure Bl. Ratio and contour plots for sources which do not require a broadened component in the FeK band. In each plot the 
data/model residuals (top panel) show the residuals which remain in the FeK band when all atomic lines have been removed and 
the reflection component is fitted with pexrav to highlight the presence of the neutral FeKa/K/3 fluorescence lines. The contour plot 
(bottom panel) show the F-test significances of the remaining residuals when the FeKcr and K/3 lines are fitted with reflionx and a 
narrow Gaussian line. The continuous outer contour corresponds to a Ax 2 = +0.5 worse fit and is intended to indicate the level of the 
continuum baseline. The closed coloured contours correspond to A\ 2 improvements of —2.3 (red), —4.61 (green), —9.21 (blue), —13.82 
(cyan) and —18.42 (magenta), which translate to F-test significances of 68%, 90%, 99%, 99.9% and 99.99%, respectively. The dashed 
vertical lines indicate the expected rest-frame energies of the FeKct, Fexxv Hea and Fexxvi Lyct transitions. Colour version available 
online. 
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Figure Bl - continued - Ratio and contour plots for sources which do not require a broadened component. 
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Figure Bl - continued - Ratio and contour plots for sources which do not require a broadened component. 
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Figure B2. Ratio and contour plots for sources which do require an additional broad component at FeK. The top and middle panels 
again show the primary residuals at FeK and their significances according to the F-test, respectively. The remaining residua after the 
addition of this broadened component (i.e., a broadened Gaussian or a diskline) arc shown in the bottom panel. The confidence contours 
and dashed vertical lines have the same meaning as in Figure B. Colour version available online. 
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APPENDIX C: MODELLING COMPLEXITIES 

There are several modelling caveats which could have an af- 
fect on the detection of FeK-band absorption lines. Here, we 
investigate whether there are any degeneracies or uncertain- 
ties in the broad-band models which could have an effect 
on the positively detected absorption lines reported in this 
work. 

CI The Soft-excess 

A total of 6/20 sources in which we have detected highly- 
ionised absorption also have a soft-excess. Thus far we have 
taken the simplistic approach of fitting any excesses below 
2 keV with a phenomenological blackbody however there are 
several alternative means of fitting this component which 
can subtly alter the broadband continuum model, and hence 
possibly influence the parameters obtained for the FeK- 
band absorption lines. To assess for this possibility we re- 
fitted the soft-excess in each of the 6 sources with two al- 
ternative models: diskbb (Mitsuda et al. f984) which inte- 
grates over the surface of accretion disc to form a multi- 
colour blackbody spectrum, and comptt (Titarchuk 1994) 
which models the excess with the Comptonised emission of 
soft seed photons by a hot corona above the accretion disc. 
Both of these alternative models have been used to good 
effect in the literature when modelling the soft-excess (e.g., 
Porquet et al. 2004; Patrick et al. 2011; Gofford et al. 2011; 
Patrick et al. 2011) 

We replaced the bbody component in each source with 
both diskbb and comptt in turn, refitted the data in each 
case, and compared the resultant continuum and absorp- 
tion line parameters to those obtained previously. For the 
comptt models a disc geometry was assumed and the in- 
put photon seed temperature was allowed to vary between 
5eV and O.lkeV, while the optical depth r and temper- 
ature T of the Comptonising plasma above the disc were 
allowed to vary freely, as was the component normalisation. 
In all cases we found that fitting the excess with diskbb 
or comptt had no significant effect on the parameters of ei- 
ther the absorption lines or the continuum, which were al- 
ways found to be consistent with those obtained with bbody 
at the 90% level. In the heavily absorbed sources such as 
NCG3783, NGC4051 and MCG -6-30-15 there were some 
small variations in column density and ionisation parameter 
for the soft band absorbers when fitting with comptt, pre- 
sumably to compensate for the fact that the Comptonised 
emission fitted with comptt can extend to harder X-ray en- 
ergies (i.e., E > 2 — 3keV) than typically seen in black-body 
derived models. Even so, because the key parameters (i.e., 
those of the continuum and the hard-band absorption lines) 
are not affected by these subtle changes we do not consider 
them particularly significant. On the whole, the parameters 
measured for the highly-ionised absorption systems in these 
sources are largely independent of the method used to pa- 
rameterise the soft-excess. 

As mentioned in Section 4.2.3 other interpretations of 
the soft-excess posit that it may be due smeared absorption 
or blurred reflection associated with an increase in optical 
depth of lowly ionised circumnuclear material. For simplicity 
we do not investigate these interpretations further as they 
likely require a restructuring of the entire broadband model 



for each source. However, we note that most the outflows 
detected in these sources have also been detected by Patrick 
et al. (2012) (see Section 6.1.3 for a detailed sample com- 
parison) who have used blurred reflection models to fit their 
Suzaku spectra which suggests that the majority of the line 
detection are robust regardless of the underlying assump- 
tions made when modelling the broad-band continuum. 

C2 Warm Absorption 

Soft band absorbers with intermediate ionisation states, i.e., 
those with log(^/ergcms _1 ) ^ 2.5 where Fe starts to be- 
come significantly ionised (Kallman & McCray 1982; Kali- 
man et al. 2004) can impart weak Fexxv and Fexxvi ab- 
sorption lines in the Fe K band which can reduce the signif- 
icance of any absorption lines associated with higher ioni- 
sation material. Absorption components with ionisation pa- 
rameters in this regime are present in 4/20 sources with de- 
tected FeK outflows. However, because all of the absorbers 
in the soft band have an implicitly assumed v out = kms -1 
throughout this work, any associated absorption lines will 
be limited solely to their expected rest-frame energies and 
only the lowest velocity absorbers in the Fe K band are likely 
to be affected. Moreover, even if this was not the case, the 
line significances through Gaussian fitting for each of the 
sources are very high, with Ax 2 /A/^ 25/2 in all cases. 
This suggests that even if there is weak Fexxv-xxvi from 
the soft-band absorbers present at higher energies, it is un- 
likely to have a significant effect on the high velocity ab- 
sorption systems detected here. For these two reasons we do 
not regard this effect to pose a problem for our detection of 
highly-ionised absorption lines in these sources. 

C3 Compton Reflection 

The ref lionx reflection model assumes a fixed face-on ge- 
ometry with all of the reflected emission originating at the 
surface of the accretion disc (Ross & Fabian 2005). This 
means that it is unable to account for the variations of the 
fluorescence line emissivity or the profile of the Compton 
reflection hump, both of which have a viewing angle de- 
pendence (Magdziarz & Zdziarski 1995), and can lead the 
model to either under- or over-predicting the Fe Kq line flux 
for a given reflection fraction if a source happens to have 
an inclined accretion disc, or has a more complex reflec- 
tion geometry. In some sources this can lead to an apparent 
hard-excess which may be a direct result of the underlying 
assumptions made with the ref lionx model rather than in- 
trinsic feature of the X-ray spectrum. Alternatively, the ob- 
served reflection could have an origin which is not associated 
with the accretion disc. For example, a strong hard-excess 
remains in MCG -2-58-22 when fitting the reflection contin- 
uum with reflionx, while Rivers et al. (2011) show that 
no such-excess is present when using the more sophisticated 
MYTorus reflection model which assumes a distant toroidal 
reflector. 

Strong hard X-ray excesses at E > 15keV have been 
reported in a number of sources in the literature in recent 
years (e.g., Turner et al. 2009; Reeves et al. 2009; Risaliti 
et al. 2009; Tatum et al. 2012). The weaker of these hard 
excesses can sometimes be accounted for by allowing the 
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iron abundance of the reflector to adopt sub-solar values, 
and thereby enhancing the prominence of the Compton- 
reflection hump relative to the flux of the FeKa fluores- 
cence line. However, this approach is unable to account for 
the observed excesses in the most extreme cases and an ad- 
ditional spectral component is required to supplement the 
neutral reflection in the hard X-ray band. There are two 
spectral interpretations which accomplish this: relativisti- 
cally enhanced reflection (e.g., Nardini et al. 2011; Reynolds 
et al. 2012) and Compton-thick partial covering (e.g., Turner 
et al. 2009; Reeves et al. 2009; Tatum et al. 2012). 

In this work we exclusively use Compton-thick partial 
covering (Nu > v^ 1 — 1-5 x 10 24 cms -1 ) to fit hard-excesses 
which cannot be accounted for through allowing the iron 
abundance of the reflionx component to vary. However, we 
note that the decision to model hard-excesses with partial 
covering is driven solely by the desire to construct a broad- 
band model that is suitable for rapidly computed Monte 
Carlo simulations, and that the presence of a Compton- 
thick partial coverer in a model does not necessarily mean 
that the source is intrinsically covered by a clumpy screen 
of Compton-thick gas. Indeed, statistically equivalent fits to 
the hard-excess can likely be achieved with models based on 
relativistically blurred reflection, but the extreme computa- 
tion time that would be necessary to run Monte Carlo sim- 
ulations on these models makes their use impractical in the 
search for highly-ionised absorption lines and we therefore 
do not consider them a viable option in this work. Even so, 
we stress that because the contribution of Compton thick 
absorption or blurred reflection is largely limited to hard 
energies, i.e., E > 10 — 15keV, the manner in which the 
hard-excess is fitted in these sources is unlikely to have an 
affect on the detection rate of highly-ionised absorption lines 
or the parameters of the associated FeK absorbers. 

C4 Partial-covering Absorption 

In addition to the Compton-thick partial covering used to 
supplement the reflection continuum in sources with hard- 
excesses, we have also used moderate columns of partial- 
covering absorption (Na < few x 10 23 cm~ 2 ) to fit spectral 
curvature and spectral variability below lOkeV. In terms of 
the sources with detected FeK outflows 12/20 have at least 
one partially covering component (including those which are 
Compton-thick), with 5 of those having more than one par- 
tial coverer in a layered or clumpy geometry. The principal 
driver behind the inclusion of partial covering absorption 
in these sources is the need to account for the variabil- 
ity in spectral shape which are apparent when more than 
one observation has been simultaneously fitted, especially 
if those observations are separated by a number of years. 
The addition of a partially covering component allowed for 
the variability observed in Mrk 766, NGC 1365, NGC 3227, 
NGC3516, NGC 3783, NGC 4051 and PDS456 to be easily 
accounted for, and enables the absorption line parameters in 
these sources to be assessed against a self-consistently fitted 
underlying continuum in all observations. It is important 
to note, however, that regardless of the underlying contin- 
uum model the detection of highly-ionised absorption lines 
in these sources is largely model independent and does not 
affect the overall statistics reported in this work. 
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APPENDIX D: FULL MODEL DETAILS 

Table Dl. Summary of continuum parameters. Notes: (1) Source name; (2) Suzaku observation ID. See table Al for details on 
stacked spectra; (3) Power-law photon index; (4) Power-law component normalisation in units of xlO~ 3 phkcv -1 cm~ 2 s~ 1 at 
IkeV; (5) Blackbody thermal temperature in units of eV; (6) Blackbody normalisation in units of xl0~ 5 X (£39/ Df ), where 
L39 is the source luminosity in units of 10 39 crgs _1 and Dm is the distance to the source in units of lOkpc; (7) Ionisation 
parameter of the ref lionx X-ray reflection component, where £ = 4nF/n; (8) Iron abundance for the X-ray reflector, in Solar 
units; (9) Reference for warm absorption parameters (if present). 'D2' and 'D3' refer to Tables D2 and D3, respectively; (10) 
Reduced \ 2 an d number of degrees of freedom (v) for the final best-fit model. 
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i ' 1 *-0.10 


1.00-'' 


D2 


1.04(978) 



* Parameter was tied during fitting. 

f Parameter was frozen at this value during fitting. 

3 Indicates that spectra were fit simultaneously. 

^ indicates that a highly-ionised outflow is detected in this source. 
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Table Dl continued - Summary of continuum parameters 



Source 


OBSID 




Power-law 
r norm 




k B T 


bbody 

norm 


Reflection 


Abs. 


x 2 M 


(i) 


(2) 




(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 2992 


stackcd[all] 


1 


oc+0.02 
BJ -0.02 


3.3512 


07 
07 


- 


- 


, c ,+0.09 
<■ Di -0.09 


1.00 f 


D2 


1.05(1610) 




stacked [be] 


1 


82* 


3.891° 


43 






2.39* 


1.00* 






NGC 3227° 


703022010 i 


1 


83 -0.01 


13.92+2 


49 
19 


- 


- 


< 1.22 


1.00 f 


D3 


1.04(5803) 




703022020^ 


1 


83* 


9.71+° 


37 
37 


- 


- 


< 1.22* 


1.00* 








703022030 3 


1 


83* 


11.4712 


37 


- 


- 


< 1.22* 


1.00* 








703022040 3 ' 


1 


83* 


4.9412 


17 


- 


- 


< 1.22* 


1.00* 








703022050^ 


1 


83* 


10.2712 


21 


- 


- 


< 1.22* 


1.00* 








703022060^ 


1 


83* 


7.6012 


21 
21 


- 


- 


< 1.22* 


1.00* 






NGC 3516° 


100031010 


2 


i i +0.03 
±J --0.03 


27.8112 


75 


- 


- 


2 '2 4 -0.10 


1.00 f 


D3 


1.01(3820) 




704062010 


2 


11* 


9.7912 


61 


- 


- 


2.24* 


1.00* 






NGC 3783° 


7010330103 


1 


QC+0.01 

BJ -0.01 


18.4812 


40 




26.40l|;g 


< 1.20 


l.()0 f 


D3 


1.04(4219) 




704063010^' 


1 


85* 


22.3912 


19 


82* 


107.72l!2;21 


1.20* 


1.00* 






NGC 4051° 


700004010 


2 


I Q + . 1 

lB -0.01 


13.2012 


07 
07 


79li 


4 33+0.56 


< 1.01 


1.00? 


D3 


1.01(4551) 




703023010 


2 


18* 


26.4312 


09 


79* 


57 Q«+ 4 - 34 


< 1.01* 


1.00* 








703023020 


2 


18* 


21.44+2 


14 
13 


79* 


33 74 + 2 -73 


< 1.01* 


1.00* 






NGC 4151° 


701034010 


2 


Q2 +0.01 
u -0.01 


22.6512 


95 
93 


- 


- 


< 1.02 


l.()0 f 


D2 


1.00(3052) 


NGC 4395° 


702001010 


1 


,.-+0.02 
DO -0.02 


1.8212 


05 


- 


- 


< 21 »ii:o2 


l.()0 f 


D2 


1.00(590) 


NGC 4593° 


702040010 


1 


64 +o.oi 


2.2012 


02 
02 


- 


- 


< 1.04 


1 07+0.19 

1D '-0.19 


D2 


1.03(1863) 


NGC 5506 


stacked [all] 


2 


r|4+0.01 


83.9812 


64 
62 


- 


- 


6.8612J 7 


1.00 / 


D3 


1.06(2376) 


NGC 5548 


stacked [all] 


1 


7-+0.01 
'°-0.01 


4.8012 


25 
25 


- 


- 


< 1.42 


1.00? 


D2 


1.05(1854) 


NGC 7213 


701029010 


1 


+ 01 
7S -0.01 


- ~~ + 
6.39^2 


02 
02 


- 


- 


~, +0 24 
5 - 01 ±0.51 


+ 47 
1.83_ 27 


- 


1.02(2686) 


NGC 7469 


703028010 


1 


oo + O.Ol 
BJ -0.01 


2.2112 


02 
02 


- 


- 


< 1.3 


1 (il+°- 16 

1D± -0.16 


D2 


1.00(2377) 


PBC J0839. 7-1214 


705007010 


1 


74 +o.oi 

' -0.01 


2.66+2 


04 






134.15+*;|? 


1.00? 




1.06(1267) 


PDS456 


701056010 


2 


47 +0.01 
^' -0.01 


3.19+2 


05 
05 










D3 


1.08(1238) 




705041010 


2 


47* 


3.07+2 


06 
06 














PG 1211 + 143 


700009010 


1 


„„+0.02 
8S -0.02 


1.1812 


02 
02 


80l 4 


1 - 1D _0.22 


< 71.87 


1.00 f 


D2 


1.04(831) 


PKS 0558-504 


Stacked [all] 


2 


30 +0.01 


13.08+2 


18 
18 


133+ 4 


10.28l2'j 7 


< 155.36 


1.00 f 


D2 


0.99(1868) 


RBS 1124 


702114010 


1 


7-+0.03 
,o -0.03 


2.7712 


06 
06 


1051- 


n Q9+ - 86 


< 27.41 


l.()0 f 


D2 


0.94(764) 


SW .12127.4+5456° 


702122010 


2 


01 + 0.02 
ul -0.02 


13.27+2 


39 
39 


2711- 


, 79 +1.01 


20ll3 


1.00? 


D2 


0.97(2063) 


TON S180 


701021010 


2 


41 + 002 
-0.02 


9.9412 


25 
25 


12412 


6 00+ - 79 
D.uu_ 77 


566.391^ + 2 


1.00 f 


D2 


1.05(922) 



© ? RAS, MNRAS 000, 1-?? 



36 Jason Gofford et al. 



Table D2. Warm absorber parameters for single-epoch sources. Notes: (1) Source 
name; (2) Suzaku observation ID; (3) Absorption zone; (4) log of absorber column 
density, in units of cm -2 ; (5) log of the ionisation parameter, in units of erg cm s _1 ; 
(6) Per cent (%) covering fraction of absorption component. All warm absorption 



components are required with Pp > 99% 


Source 


OBSID 


Zone 


logJV H 


log? 


fcov 




C91 
l+J 


(J) 


I 4 ; 






1H 0419-577 


stacked [all] 


1 


24.33l°;°! 


u - o± -0.08 


65l| 






2 


23.27±g : gi 


9 n o+0.01 


17l 2 


3C 382 


702125010 


1 


91 m+ 003 


9 70+0. 16 
z - '•-'-0.31 


100' f 


3C445 


702056010 


1 


90 90+0.01 

zo.zo_ 01 


1 19+0.07 
1 - lz -0.16 


£ 98 


4C +74.26° 


702057010 


1 


21.571°;!° 


1 61 + - 08 


100-^ 






2 


20.96t°° 5 5 


neutral 


100-f 


APM 08279+5255* 


stacked [all] 


1 


99 R7+ ' 09 


neutral 


wo f 


Ark 564 


702117010 


1 


o t i ^-4-0 01 

24 - 48 ±o.oi 


o i ^-t-0 04 

9 4fi^ 
Z -* D -0.19 


41 + 1 






2 


22.731°,;°,?, 


05+010 
u - yJ -0.09 


lOO-'' 


CBS 126° 


705042010 


1 


22.48l°;°| 


9 Q-+0.02 
-0.02 


100 / 


ESO 103-G035° 


703031010 


1 


-i-d 09 
23.07t°'° 2 


neutral 


98+1 






2 


_i_n no 

23.02t°;° 2 


-i_n oq 

0.82l°;i 


98* 






3 


91 74+0. U4 
z± -'^-0.04 


neutral 


lOO-f 


IC 4329A 


stacked [all] 


1 


, in 1 7 

91 17+ ui ' 
z± - i ' -0.09 


l.98l°;! 


100 f 






2 


21 94+ 001 
-0.01 


9 7 +0.01 
1 -' I '-0.01 


ioof 






3 


91 vr;+0.05 
zi - '°-0.05 


, 24 +0.13 


100 f 






4 


91 1 Ej+0.03 


9 cr9+0.07 
z -° -0.03 




IGRJ21247+5058 


702027010 


1 


99 54+° 13 


neutral 


< 21 






2 


21.94l°;° 2 


neutral 


ioof 


MR 2251-178° 


704055010 


1 


20.98t°;J° 9 




100 f 






2 


21-43±°:fg 


9 91 +0-08 
z - zl -0.08 


100 f 


Mrk79 


702044010 


1 


23.19«° 3 3 


9 -1 o+0.07 
z - ±o -0.09 


37 + 6 






2 


91 R9+ 004 


1 44+0- 06 
-"-• 44 -0.04 


100 f 


Mrk205 


705062010 


1 


> 24.68 


neutral 


59 + 2 


Mrk279° 


704031010 


1 


20.68l°;i 2 


1 26+ ' 28 


100 f 


Mrk335 


701031010 


1 


ZU.^O_Q Q2 


„ ^o-t-0 20 

0.52l°; 2 ° 


100 f 


Mrk841 


stacked [all] 


1 


20.39t°; 3 24 


1 84+ 041 


100 f 


MCG -2-58-22 


704032010 


1 


_ , , _ _i_o ,; i 
Z4t.lO_ 04 


, in on 

9 O Q \ U. iV) 

Z -JJ_0.13 


35 + 2 


MCG -5-23-16 


700002010 


1 


24.181°:°! 


1 71 +0.06 
1- ' -0.06 


49+1 






2 


22.681°;°! 


0.58l°;°! 


lOO-f 


MCG -6-30-15° 


stacked [all] 


1 


90 70+0.02 
zo.(o_ 02 


9 91 

z * Zi -0.05 


< 18 






2 


91 49+O.02 


n 79+0.02 
u * ' z -0.03 


100 f 






3 


^ ~--.+0 01 

91 7Q^ U,U 
zl - ly -0.01 


_ .-l-O 02 
A -'^-0.02 


100 f 


MCG +8-11-11 


702112010 


1 


20 59+ - 05 
zu.oz_ 02 


< 0.18 


100 f 






2 


91 02 + °' 13 


9 q c;+0. 10 
Z - JO -0.14 


ioof 


NGC 1365° 


705031020 


1 


. , . , LA 1 O 

24.17l°; 3 
23.791°;°! 


neutral 

1 09+0. 06 
J - oz — 0.07 


31 + 1 
> 0.96 


wrr 9qQ9 




I 


91 S0+ ' 02 

zi.su_ n 02 
22.1241°;^ 


Ulrl Ojl 

1.79 + 0.08 


95 ± 1 
95 f 


NGC 4151° 


701034010 


1 


99 91 +0.06 


z -°'-0.18 


^ 98 






2 


23.001°°! 


fis+ 02 
U.DJ_ 02 


> 98* 






3 


23.691°;°! 


1 8fi+°- 07 
1 - Sb -0.06 


68 + 5 


NGC 4395° 


702001010 


1 


23.9ll°;° 8 6 


neutral 


38 + 2 






2 


22.34l°;° 3 


56+ 014 
U - OD -0.13 


83+1 






3 


21.58l°;! 7 


1 64+°- 17 
J -° -0.27 


lOO-'' 


NGC 4593 


702040010 


1 


21.99l°; 2 | 


2 96+°- 18 


100^ 






2 


21 70+0.24 


9 90+0. 06 
^•^8-0.09 


100 / 



* indicates a parameter was tied during spectral fitting; 

f Indicates a parameter was frozen during spectral fitting; 

° indicates that a highly-ionised outflow is detected in this source. 
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Table D2 - continued - Warm absorber parameters for single-epoch sources. 



Source 


OBSID 


Zone 


loff Ntj 
iu & 1 " H 


loef 


f 

J cov 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 5506* 


stacked [all] 


1 


21.51+2-22 

—0.02 


neutral 


100 / 






9 
Z 


99 QQ+0-03 
zz.oy_ 02 


i 

± " oo -0.07 


, QQ 






3 


22 29+2'2o 
— 0.02 


neutral 


> 99 






4 


24 m +0 01 


11(_- U.L1 ell 


35 i 1 


NGC 5548 


Kt p\ pWoH fa 1 1 1 


1 


22 15+°-?I 

ZZ. 1J_q 12 


2 78+2 22 


100-f 






2 


21 44+2' 2 ,i 
-0.05 


2 20+2'?? 


> 92 






3 


21 84+ ' 03 
zl -°^-0.03 


38+ 12 


46 ± 2 




70901 ^010 


i 


91 K9+006 


neutral 


i nn/ 

1UU 






2 


21 94+0-03 
' i± ' t ' -0.03 


i 99+0.15 

J--zz_ q 12 


100 / 


NGC 7469 


703028010 


1 


91 Q9+010 
z±.oz_ og 


9 ocr+0.07 
z.oj_q 08 


100 / 


PG 1211+143 


700009010 


1 


22 4Q+0- 13 


9 SS+0- 07 
^•»»_0.09 


100^ 


PKS 0558-504 


stacked [all] 


1 


24.28+ ,; 04 


1 O4+0-09 


38 ± 2 


RBS1124 


702114010 


1 


24 62+0' 22 


neutral 


56 ± 1 


TONS180 


701021010 


1 


24 14+0 03 


1.92+°,;,° 


64 ± 1 


SWJ2127.4+5456<> 


702122010 


1 


91 i 9 +0.03 
z± - iz -0.03 


0.83+° 5 5 ° 


100-f 
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Table D3. Warm absorber parameters for jointly fit spectra. Notes: (1) Source name; (2) Suzaku observation ID; (3a, 4a, 5a, 6a) logarithm of absorber column density, in units 
of cm -2 ; (3b, 4b, 5b, 6b) logarithm of absorber ionisation parameter, in units of ergems -1 ; (3c, 4c, 5c, 6c) Absorber covering fraction, in per cent. All absorption components are 
required at Pp > 99% significance. 



Source 


OBSID 


logJV H 


Zone 1 
log? 


/cov 


log^H 


Zone 2 
log? 


/cov 


logAf H 


Zone 3 
log? 


/cov 


logJV H 


Zone 4 

log ? /cov 


(1) 


(2) 


(3a) 


(3b) 


(3c) 


(4a) 


(4b) 


(4c) 


(5a) 


(5b) 


(5c) 


(6a) 


(6b) (6c) 


3C111 


703034010 
705040010 
705040020 
705040030 


21.83±g;°i 
21.83* 
21.83* 
21.83* 


neutral 
neutral* 
neutral* 
neutral* 


wo f 

100* 
100* 
100* 


















3C120 


700001010 
stacked [bed] 


20.961°,;°,! 
20.96* 


neutral 
neutral* 


100 / 
100* 


















Mrk 509 


705025010 
stacked [be] 
705025010 


21 40+°'° 4 

21.40* 

21.40* 


2 1Q+ - 04 
z - ly -0.05 

2.19* 
2.19* 


wo f 

100* 
100* 


















Mrk 766° 


701035010 


■"•OO-0.02 


z - ±a -0.03 


94 + 4 


21.86t°;° 6 6 


i sn+0' 04 


100-f 














701035020 


22.66* 


2.19* 


< 14 


21.86* 


1.80* 


100* 












NGC 1365° 


702047010 


24 1 1 1 +0.02 
Z ^-- L1 -0.03 


9 49+ - 24 

-0.24 


55+5 


23 02+ - 01 


neutral 


42± 4 


99 sq+0.06 
zz.oz_ 06 


o 91 +0.04 
°- z -0.04 


42* 








705031010 


24.11* 


2.42* 


40+J 


r>o CO + 0.01 

zo.oo_ 01 


neutral* 


58 ±2 


23.20l°;°| 


3.21* 


58* 






NGC 3227° 


703022010 


, „^o m 

99 1 fi+ u - ui 


0.81 ± 0.05 


W0 f 


23-06l°;°i 


2.34 + 0.02 


55 ± 5 














703022020 


22.16* 


1 oo+°-° 5 


100* 


23.06* 


-0.11 ± 0.04 


86 ±8 














703022030 


22.16* 


J-Ci 1 1 

67 + 
U - D '-0.11 


100* 


23.06* 


1.53 ±0.05 


87 ±4 














703022040 


22.16* 


-, 1O+0.09 
1 - ±o -0.08 


100* 


23.06* 


1 ie+0.19 
1 - ±o -0.18 


83 ±6 














703022050 


22.16* 


6fi+ 016 
U - DD -0.16 


100* 


23.06* 


1.24 ±0.04 


91 ±4 














703022060 


22.16* 




100* 


23.06* 


o.67t;; ; i? 


90 ±5 












NGC 3516° 


100031010 


94 22+°'° 2 
z4 - zz -0.02 


2 00+ 02 


46 + 4 


23.27+o. i 


1.89 ±0.04 


> 93 


91 „ 9 +o.oi 


1 01+ 007 


> 99 


21 50+ 010 


< 0.14 > 99* 




704062010 


24.22* 


2.00* 


48l 9 8 


23.27* 


1.89 ±0.04 


58 ±11 


21.32* 


1.01* 


> 77 


21.50* 


< 0.14* > 77* 


InOC 3 (oo v 


7n 1 HQ QO.1 c\ 

(U1033010 


9q 7R+0.02 
23.76I 02 


i Qo+0.04 
1-y -0.12 


27 i 3 


99 9 o+0.01 
zz.zo_ 01 


1 r;K+0.01 
L56 -0.01 


100^ 


99 in+ 09 
ZZilu -0.09 


9 77 +0.07 
z - 1 1 -0.07 


100^ 








704063010 


23.76* 


2 4i+ - 06 


22 + 1 


22.23* 


1.92* 


100* 


22.10* 


2.77* 


100* 






NGC 4051* 


100004010 


24 S6+ ' 01 


r, 1 1+0.02 

z ' 1J --0.02 


26 + 2 


2 o u +0.02 


1 65+ ' 05 


< 14 


21 41+ - 01 

^■^l-O.Ol 


1 4R+003 


100? 








703023010 


24.36* 


2.11* 


50 + 6 


23.14* 


1.65* 


40 ±2 


90 54+0-05 


1.46* 


100* 








703023020 


24.36* 


2.11* 


< 10 


23.14* 


1.65* 


32 ± 3 


21.41* 


1.46* 


100* 






PDS456* 


701056010 


22.98t°;°3 


n 1 rr+0.04 
z - ±J -0.04 


28 + 2 


22 00+ ' 08 


O 1O + 0.10 

0,10 -o.ii 


100-f 














705041010 


22.98* 


2.15* 


56 + 3 


22.00* 


1 65+ ' 11 


100* 













* indicates a parameter was tied during spectral fitting; 
f Indicates a parameter was frozen during spectral fitting; 



indicates that a highly-ionised outflow is detected in this source. 
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Table D4. Principal omission components in the 5-10 keV band. Notes: (1) Source name; (2) Suzaku observation ID; (3) Emission 
line identification. 'Broad(GA)' and 'Broad(DL)' indicate that a broadened FeK emission component, which was fitted with 
either a broadened Gaussian or a disk-line, respectively; (4) Measured line energy in the source rest-frame, in units of keV; (5) 
Measured line energy-width, in units of eV. Unresolved lines were fit with widths fixed at 10 eV; (6) Line intensity, in units of 
X 10 — 6 photons cm~ 2 s —1 ; (7) Line equivalent width, in units of eV. For the FeKo line the EW of the corresponding FcK/3 line 
is given in square brackets ; (8) Change in fit statistic (and degrees of freedom) when line is removed from the best-fit model; (9) 
Line significance according to the F-test. 



i rrf 

UU HI l^>_> 


OBSID 


Line TD 


Pri r>t'o'"\," 


\J W 1U LI1 




EW 


Ay 2 / A,/ 




F 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


( 7 ) 


(8) 




(9) 


ill u4iy-o 1 1 


stacked [all] 


rC 1\G. [IvpJ 


f. QQ+0.02 


in/ 


4 - b -l.l 


01 +4 r** j- 1 1 


oO. 4 


> 


no nn 

yy.yy 


Q f~i 1 1 1 
Ob 111 v 


/ UoUo4UlU 


re l\o.|l\_pj 


fi , Q +0.02 
u - oy -0.02 


in/ 


17 1+ 3 ' 1 


0-7+8 fio~l"2i 


1 no 9 /9 
lUZ.z/ z 


> 


nn nn 

yy.yy 




/U0U4UU1U 


XT',, IV" ,^ OV /Ql 

rC lVQ![lvpJ 


R OO+0.02 


in/ 
1U J 


Ofi 1+4.0 
26-1-4.0 


4b_ 7 [0 ± 1J 


93.0/2 


> 


99.99 






Fc xxv I ^ 


fi S5+ ' 10 


10^ 


°-°-3.8 


L( -6 


13.0 




99.99 




/ U0U4UUzU 


rC !VQ![lvpJ 


fi ,0+0 04 


in/ 


14 o+4.4 
14.3_ 4 4 


19_ 6 [4± lj 


/i n 

4y.u 


> 


nn nn 

yy.yy 




7Arri / | noun 
/UoU4UUaU 


TJV, TV OV (Ql 

tc Ka|KpJ 


R -n+0.03 
6 - 40 -0.03 


inf 


on k+ 4 - 4 
20.b_ 4 4 


29_ 6 [5 ± 1J 


68.7 


> 


99.99 


q pi i on 
DO IzU 


7nnnni 01 o 

i UUUU1U1U 


rC lVQ![lvpj 


n OQ + 0.03 

o-^-o.oa 


in/ 


is n+ 4 - 4 
lo.U_ 4 4 


o,+8 rc + li 

31_ g [£>_iJ 


^7 ^ /9 

o/.O/z 


> 


no 00 
yy.yy 








fi K/l+0.12 
b - 54 -0.12 


o nn +162 
300_ 86 


or) + 8.8 

-8.8 


co+36 
O»_ 40 


l0.4/ O 




no fi9 
yy .oz 




stackcd[bcd] 
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indicates that a highly-ionised outflow is detected in this source. 
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Table D4 - continued - Principal emission components in the 5-10 keV band 
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stacked [all] 


FeKa[K/3] 


6.38i° ° ; 


r,+15 

0i -16 


23.412.? 


noli" [iTtl] 


439.5/3 


> 


99.99 


NGC 7213 


701029010 


FeKo[K/3] 


6 4n+ - 01 


< 38 


20.3±11 


79l 9 9 [1211] 


270.0/3 


> 


99.99 






Fexxv 


6 64+ ' 04 
°' D4 -0.03 


lO-f 


5 2+ 1 ' 8 


19l 6 6 


22.5/2 


> 


99.99 






Fc xxvi 


6 95+ - 03 
°' ao -0.03 


10/ 


7 - 5 -1.8 


2^+8 


43.2/2 




99.99 
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The Suzaku view of highly-ionised outflows in AGN. 



Table D4 - continued - Principal emission components in the 5-10 keV band 



Source 
(1) 


OBSID 

(2) 


Line ID 

(3) 


Energy 
(4) 


cx-width 
(5) 


Intensity 
(6) 


EW 

(7) 


A x 2 /Au 
(8) 


Pf 

(9) 


NGC 7469 

PDS456 
PG 1211+143 

SWJ2127.4+5654^ 


703028010 

701056010 
700009010 

702122010 


FcKa[K/3] 

broad(GA) 

blend? 

FcKa[K/3] 

Fexxv 

FcKa[K/3] 

broad(GA) 


6 40+ - 01 

6 41+ 012 

7Q4 +0.17 
' - U4 -0.15 
g o 7 +0.04 

6 64+ - 03 
6 40+ - 07 
6 22+ - 14 


< 65 

262+ 142 
zo -109 

QOO + 214 

zoo -126 
qo+46 


23.1+3 2 
20.5+e* 

— o.o 

3-211:1 
3-211:1 
6-911:? 

6 0+ 4 ' 3 

39 1+ 12 ' 9 
— 12.1 


941:11 [911] 

84+26 
°*-24 

72l| 2 

59l 2 ° [10l 3 ] 
130+ 34 . 
16111 PHI 
lOOlf, 3 , 


253.4/3 
23.1/3 
12.8/3 
20.2/2 
50.7/3 
10.8/2 
23.6/3 


> 99.99 

> 99.99 

> 99.22 

> 99.99 

> 99.99 

> 99.99 

> 99.99 
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